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PREDICTION OF UNSTEADY AERODYNAMIC LOADINGS CAUSED BY LEADING 
EDGE AND TRAILING EDGE CONTROL SURFACE MOTIONS IN SUBSONIC 
COMPRESSIBLE FLOW— ANALYSIS AND RESULTS 

By 

W. S. ROWE, M. C. REDMAN 
F. E. EHLERS and J. D. SEBASTIAN 
Boeing Commercial Airplane Company 

SUMMARY 

A theoretical analysis and computer program has been developed for the 
prediction of unsteady lifting surface loadings caused by motions of leading 
edge and trailing edge control surfaces having sealed gaps. The final form 
of the downwash integral equation has been formulated by isolating the singu- 
larities from the non-singular terms and using a preferred solution process 
to remove and evaluate the downwash discontinuities in a systematic manner. 
Comparisons of theoretical and experimental pressure data are made for several 
control surface configurations. The comparisons indicate that reasonably 
accurate theoretical pressure distributions and generalized forces may be ob- 
tained for a wide variety of control surface configurations. 

INTRODUCTION 

Historically, much effort has been expended to adapt conventional oscil- 
lating lifting-surface solution techniques to wings with control surfaces. 
Extensive use has been made of the reverse flow theorem (reference 1) to con- 
struct an analytical upwash function that is used in reverse flow to produce 
"equivalent" generalized forces. The generated upwash function is obtained 
by the expedient of smoothing over, in one way or another, the slope discon- 
tinuities and corresponding singularities within the pressure functions. The 
"equivalent" upwash function does not match the control surface slope accur- 
ately at all points, but does produce reasonable "indirect" generalized forces. 
The usefulness of this method deteriorates when the "direct" hinge-moment term 
is to be determined. 


The work of reference 2 is one of the few subsonic methods developed 
to determine the 'direct' surface loadings using pressure terms that are 
capable of correctly representing the known singularity functions around 
the boundaries of the control surface. However, it was found that solu- 
tions obtained from the method of reference 2 were highly sensitive to the 
relative location and number of control point collocation stations used in 
the analysis. The sensitivity may be attributed to the particular solution 
process being applied that assumes that discontinuous downwash distributions 
may be approximated by a linear combination of polynomial represented down- 
wash sheets that satisfy the boundary conditions at a select set of control 
points.' Changing the control point locations by relatively small amounts 
results in large changes in the unsteady 1 oadings. consequently , the method 
requires calculation of downwashes at many stations and seeking solutions 
in a least-squares-error sense* 

Solution sensitivity remained until pressure distributions were de- 
veloped (reference 3) that identify the functional distribution and singu- 
larity strengths required to produce identical mathematical downwash dis- 
continuities contained in the kinematic distribution. A preferred solu- 
tion process was developed by subtracting the discontinuous mathematical 
downwash distribution from the discontinuous kinematic distribution re- 
sulting in smooth downwash distribution for which standard lifting-surface 
solutions could be applied. The resulting loadings would then be relatively 
insensitive to .the locations and number of control points used in the 
analysis. 

The present work represents an extension of the analytical methods 
suggested in reference 3 to provide a capability of numerically predicting 
thesinsteady loadings caused by control surface motions that is relatively 
insensitive to locations of collocation stations on the surface. Documen- 
tation of the computer program design, usage, and limitations is provided 
in reference 4. 



SYMBOLS 


All quantities are dimensionless except as indicated. 


a nm 

b [length] 
b Q [length] 
c 

A 

c 

f(n) 

9(5.n) 

Hj(x,y) 

i 

K(x,?,y,n) 

k 

k 

M 

M 

N 

n 

P [Force/area] 
P u [Force/area] 
AP [Force/area] 


Unknown coefficient of assumed pressure modes, equation 69 
Local semi chord 
Reference length 

Local chord length nondimensional with respect to bg 
Variable used in spanwise pressure modification function, 
equation 49 

Spanwise pressure difference distribution, equations 56 
and 69 

Chordwise pressure difference distribution, equations 56 
and 65 

Mode shape of mode j (see equation 5 for dimensions) 

Kernel function, equation 6 

ubg 

Reduced frequency = -y- 
Index variable, equation 73 
Mach number 

Number of downwash stations on a downwash chord 

Number of downwash chords on semi span 

Index variable, equation 69 

Pressure on lower surface 

Pressure on upper surface 

Pressure difference, P.-P (positive upward) 

Xr U 


<1* 


AP. [(Force/area) /unit q.] Surface pressure difference in mode j 

J J J 


Q. . 

W 1J 

Q?. 

V 1J 


R 
S 
s 
V 

¥ 

*3 

x,y,z 

_x 

x m 

m 


[length] 
[length/ time] 


cs 


is 


Generalized coordinate amplitude for mode j 
Generalized forces (see equation 77 for dimensions) 
Sectional generalized forces (see equation 76 for 
dimensions) 

W + B V 

Semi span 

Nondimensional semi span, S/b Q 
Free stream velocity 


Kinematic angle of attack or nondimensional normalwash 

w 

V 

i cot 


"/'’j 6 

Coordinates nondimensional with respect to b Q 
Local chordwise coordinate (x - x m )/(b/bg) 
Nondimensional coordinate of the mid-chord line 
Nondimensional coordinate of the leading edge 
Nondimensional coordinate of the trailing edge 
x - 5 

Streamwise coordinate of hingeline at side edge of 
control surface 

Streamwise coordinate of planform leading edge at side 
edge of control surface 

Local coordinate for full chord control surface, 
equation 46 

Local coordinate for trailing edge control surface, 
equation 42 



y 0 

y s 

z s 

B 

e h 

AC P 

A c [radian] 

[radian] 

5 c 

5»n»C 

© H (n) [radian] 

x 

p [mass/length 3 ] 
t [time] 

w [l/time] 

0 

* 


y-n 

/ 

Spanwise coordinate of control surface side edge 
Surface deflection (see equation 4 for dimensions) 

Vl-M 2 

+ tan 2 A c 
■^B 2 + tan 2 A fc 

Non-dimensional pressure difference between lower and 
upper surfaces (P^ - P u )/( H pV 2 ) 

Sweep angle of the control surface hingeline, positive 
swept back 

Sweep angle of leading edge, positive swept back on right 
hand side 

Hingeline coordinate at span station n 
Dummy variables for (x,y,z) 

Rotation angle of control surface hingeline at (n)> 

measured in the plane perpendicular to the y-axis and 
positive trailing edge up 
Multiplication symbol 
Density of the fluid 
Time 

Circular frequency of oscillation 
Chordwise angular coordinate, equation 69 
Spanwise angular coordinate, equation 69 
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ANALYTICAL FORMULATION 
The Integral Equation and Its Elements 

This report describes the procedures used in calculating generalized 
aerodynamic forces existing on a lifting surface caused by motions of leading 
edge or trailing edge control surfaces in subsonic flow. The analysis is 
applicable to planform configurations having full span or multiple partial 
span control surfaces located anywhere along the leading or trailing edge of 
the planform. 

The analysis coordinate system is defined in figure 1 for a typical 
trailing edge and/or leading edge control surface configuration: 



Figure 1 - Coordinate System Definition 
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Solutions are obtained by evaluating the downwash-pressure Integral 
equation that has been derived from the linearized potential equation (1). 

v 2 *(x,y,z,t) - — [v ^ — + — 1 ♦(x.y.z.t) - 0 (1) 

a 2 L ax ot J 

The downwash- pres sure Integral equation derived In reference 5 provides 
the basis for the analytical formulations used In developing numerical methods 
for predicting unsteady loadings on lifting surfaces. 


The general form of the integral equation derived In reference 5 and 
used throughout this discussion is given by the expression 
w„ i 


4wpV 2 


j j APU.n) K(x.e.y.n) d(b Q e) d(b Q n) 


( 2 ) 


where 


V 


represents the kinematic downwash caused by motions of the lifting 
surface 


AP(5,n) represents the pressure difference between the upper and lower 
surface and is the only unknown in the equation 

K(x,5,y,n) represents the kernel function that defines the normalwash 
at x,y due to a pulsating pressure doublet at £.n. 

Solutions of the integral equation are obtained by assuming that the 
lifting surface deflection shape, Z s , may be represented by a finite sum of 
discrete mode shapes that are deflecting In a sinusoidal manner. That is, 
the deflection shape is assumed to be represented by 

h - E V x ' y) q j e '“ t (3 

j 

where Z $ has dimensions of length and is a function of time. 

H.(x,y) is the deflection shape of the j-tfi mode 

J 

is the generalized coordinate defining how much of the jth mode is 
contained in Z $ , 





The dimension of the product Hj(x,y)qj Is In units of length. That Is, 

If Hj(x,y) has a dimension of length, then Is non-dimensional, or If 
Hj(x,y) Is non-dimensional, then q^ has units of length. 

The kinematic downwash term of equation (2) Is developed from the surface 
motion definition and Is given as 


DZs 

Dt 


+ vi— 1 
l »t »(b 0 x) J 


*■?'[ 


aHj(x.y) 

a(b n x) 




■]' 


.1«ot 



(4) 


The total pressure difference between upper and lower surfaces is also 
dependent upon the modal displacement and is given as 


AP - P £ - 


P 


u 


» ■ 2>j) ij e1 “‘ 
j 


(4a) 


where AP. 

J 


3(AP) 

a q. 


(1 .e. , aP, is per unit q .) 
J J 


The final form of the downwash Integral equation is obtained by 
Inserting the expressions for the kinematic downwash and pressure difference 
Into equation (2) and removing the qje^ ut term from both sides of the 
equation that results in the expression 


3 H. icoHj 

3(b Q x) + V 


4lTpV Au/t-face 


J J AP.-U,n) K(x,5,y,n) d(b Q e) d(b Q n) 


(5) 
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The pressure distribution APjU.n) Is the only unknown function within 
the Integral equation and Its dimensionality Is established by noting that 

the product of K(x,5,y,n) d^Od^n) Is non-dimensional and (l/4irpV 2 ) has 

dimensional units of ( l/( Force/area ) ) which requires that the dimensions of 

AP - (^,n) is equal to the product of (Force/Area) times the dimensions of the 

kinematic downwash term. 

That Is, If H.j(x,y) Is dimensional, then the kinematic downwash is non- 
dimensional and results In having APjU,n) take on dimensions of (Force/Area). 
However, If Hj(x,y) Is non-dimensional, then the kinematic downwash has 
dimensions of (1/length) which results In having dimensions of APjU.n) 
defined as (Force/Area) /length . 

The subscript (j) Is omitted In the remaining sections, up to the 
section on generalized forces, with the understanding that the pressure terms 
being determined are those terms necessary to satisfy the jth mode shape 
boundary conditions. 


The K(x-5,y-n,M,k) term of equation (2) represents the kernel function 
of the integral equation that describes the surface normalwash at (x,y) due 
to a pulsating pressure doublet located at (£,n). The kernel function is 
defined as 


K(x-S,y-n,M,k) = e 


-if(x-C) 

= A 1 


lim 

z+0 


9 2 

9T 


x-£ 


/ 


ex Pvfr ( 

^x-M[x 2 +e 2 (y- n ) 2 +0 2 z 2 _ 

\) 


[* 2 +e 2 (y-n) 2 +B 2 z 2 ] 



dX 

( 6 ) 


Reduction of equation (6) into a form easily evaluated by routine 
numerical procedures has been accomplished by Watkins, Runyan, and Woolston 
(reference 6) and presented in non-dimensional form by considering the vari- 
able as being referenced to some chosen length (b Q ) and introducing the 
reduced frequency parameter k * b Q u/V the function becomes 


9 


/ 


K(x- 5 ,y-n.M,k) = b Q2 y Q ° [ k|y 0 l [ K l |ky 0 ,+i i' ( J 1 l ky ()l" L l I ky o0 ' 1 

h -i 

-i J*[ t /(1+t 2 )^] exp (1k|y 0 Jx) dr 

- x 0 /(x 0 ! +6 2 y 0 2 ) % -exp(i|- [x 0 -H(x 0 2 +e 2 y 0 2 ) % ])j 

h " [v M(x o !+62y o 2)% ]/ B2y o 


(7) 


where 


Singularities contained in equation (7) have been identified and are 
isolated in the following form 

K(x»£»y»ri) s k ns (x *y Kg(x»5,y,n) (Q) 

where K ns (x,C»y*n) is the non-singular part of the kernel function. 


The singular part is identified as 


e “ ikx Q f 

K s (x,cy,n) ■ -gT— yjr [ 


Hx„/(x„ 2 +6 2 y„ 2 ) % ] + 1k/(x 0 2 ts 2 y 0 2 )' 


0 ' ' 0 ■’o 


- |l £„|(x 0 2 *6 2 y 0 2 ) % -* 0 l - gr [x 0 /(x 0 2 +e 2y 0 2 ) % - m] j 


(9) 


This expression contains the dipole, inverse square root, and logarithmic 
singularities identified in reference (6). The last term of equation (g) is 
not singular. However, it is included since special means are required in 
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j 


evaluating the downwashes in regions near the downwash chord that are produced 
by this term. 

Insertion of the above definitions Into the downwash integral equation 
results in the following expression that is to be evaluated over the surface 
to satisfy the boundary condtions. 


4irpV 2 If = ff AP(5.n) K ns (x,£y,n) d(b 0 Od(b 0 n) 

-jfjfAPU.n) e' ikx o ^ [l+x 0 /(x 0 2 +e 2 y 0 2 ) % ]dcdn 
+ Ik ff AP(C,n) e y(* 0 ! +6 2 y 0 ! ) % dWn 

if 2 ff -ikx . 

- £ JJ AP U,n) e 0 -en|(x 0 2 +B 2 y 0 2 r* - x Q | d£dn 

‘ wff APU,T]) e lkX ° [ x o /( V +e2 V )JS ‘ M ] d5dn (10) 


The first integral represents the contribution to the downwash at (x,y) 
due to the non-singular part of the kernel function. The general shape of 
KpgU.f^y.n) is usually quite smooth for low Mach numbers and is due to the 
strong sinusoidal characteristic of e ikx 0 in the expression. However, the 
chordwise characteristics change quite rapidly for Mach numbers greater than 
M = 0.85 as is demonstrated in figure 2 that represents the chordwise variation 
in the imaginary part of K ns (x,f;,y,n) with Mach number. The (x-£) coordi- 
nate of figure 2 represents the chordwise difference between the pulse sending 
point located at (C,n) and the downwash station (x,y). 
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Since the functional distributions of K n$ are quite different forward 
and aft of the downwash station the chordwise integration Interval is sub- 
divided into two intervals x^<£<x and x<£sx t to provide an accurate 
evaluation of the chordwise integral. 



Figure 2 - Chordwise Variation of K ^ as a Function of 
Mach Number for y-n=0.05, k=0.7854 


The second, third, and fourth integrals of equations represent the 
downwash contributions due to the dipole, inverse square root, and 
logarithmic singularity terms respectively. The last integral does not 
contain a singularity in the spanwise direction, however, special procedures 
are required to evaluate its downwash contribution since a finite 
discontinuity exist for values of n-^y* 
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Evaluation of the Dipole Term 


The downwash contribution due to the dipole term is evaluated using the 
method suggested by Hsu (reference 7). The dipole integral of equation (10) 
is given as 

^ ^t ™ i kx r t 

fr p =-// 4P(C,n)e 0 [l+(x-{)/((x-{) 2 ^ (y . n) 2)%J d5dn ()1) 

-s x, 

The downwash contribution due to the dipole term is determined by 
considering the normalwash at a point (x,y,e) that is located just off 
the surface shown in figure 3. 



Figure 3 - Coordinate System Used in the Kernel Function Derivation 


The kernel function describing the normalwash at (x,y,e) due to a 
pressure pulse at (S.n.O) is given by 


K(x,£,y,n) = e 




( 12 ) 
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where 


j 


l 


h 

E 0 = f e ikX (r 2 +X 2 )" % dX 

-oo 

h = (x 0 -MR)/3 2 ; R = (x 0 2 +B 2 r 2 ) % ; x Q = x-£ 


The dipole portion of the kernel is 
. i kx 

K dp = + e 0 (1 +Xq/R) (-1 + 2 sin 2 ip) 

Using the geometric properties of figure 3 that define 
r 2 = r Q 2 +e 2 ; r Q = y Q = y-n; sin ^ R = (x Q 2 +3 2 r 0 2 )^ 

The dipole portion then takes on the definition 



< r o 2 - e ‘) 

( r Q 2 H *) 2 


(i + x 0 /r 0+ o( £ >>) 


03 ) 


where the terms represented by 0(e 2 )-*0 as e-K). 

The normalwash at (x.y.o) is obtained by a limiting process allowing 
e -H) after the normalwash integral has been evaluated as indicated by the 
following: 


s 

W dp (x,y,0) = lim - J 
e+0 -s 


(y 0 2 +e 2 ) 



-4P_(C»n) 

(s 2 -n 2 )^ 



[l+x 0 /R 0 ]d?dn 


(14) 


Following the development of reference 7, let G(x,y,n) be the chord- 
wise integral 


14 


(Hx 0 /R q ) d? 


(15) 


S(x,y,n) - J* e' ,kX ° 

J (s 2 -n )"* 


and as my this integral becomes 


lim G(x,y,n) ■ G(x,y,y) 

my 


One of the singularities is then subtracted from equation (14) that 
results in the expression 


H dp (x,y.o) - 11m f y 0 dn 

e. 0 l ' ° (y 0 + E J ) y n 


1 (v Z -£*) 

G(x,y,y) f — 6 (s 2 -n 2 )^ dn 

_ s (y 0 2+ e a ) 


) 


(16) 


The finite part of the second integral is 


lim - / 


(y 2 +e 2 ) 
e+o -s vjr Q ' 


- y ° — -- (s 2 -n 2 ) % dn - n 


The first integral in equation (16) is also singular. However, it may 
be evaluated using Cauchy integrals provided that the term 


[G(x,y,n) - G(x,y ,y)] 
y - n 

exists in the limit as n+y and takes on a definition of the first derivative 
of G(x,y,n) with respect to n* 
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The need of applying Cauchy integral procedures to equation (16) may be 
circumvented by subtracting an additional singularity. 

Defining 

H(x,y,n) « (s 2 -n 2 ) [G(x>y - ,r] j - I (17) 

and let 

H(x,y,y) = lim H(x,y,n) 

W 

» - (s! - y2) I? ts(xof.Ti)O n . y 


then 


W dp (x,y,0) = lim 
e->0 



■-V-1— I LH(x>y»o) - H(x t y,y)3 d 

°(y 0 2 +e 2 ) 2 ( s*-n*) % y ‘ n 


- H(x,y,y) 



(y 0 2 -e 2 ) 1 1 

_Ji drt l + 

(y Q 2 +e 2 ) (s 2 -n 2 )^ J 


TrG(x,y,y) 


08) 


The second integral limits out to zero as e-*-0 and equation (18) takes 
on the integral form of 


s 

W^ P (x,y,0) = - f 
-s 



LG(x.y,n)-G(x,y > y)] 

(y-n) 2 


(s 2 -y 2 )G^(x,y,y) . 

y-n 1 


+ Tr G(x,y ,y) (19) 

It should be noted that, as e+y, the integrand of equation (19) tends 
to take on the definition of a second derivative of G(x,y,n) with respect 
to n» 



However, the second derivative does not exist since the integrand is singular 
at the downwash chord n = y. An example of this singular characteristic when 
n-*-y is shown in figure 4 which represents a plot of the spanwise integrand 
of equation (19) applied to the analysis planform shown in the sketch. 



Figure 4 — Spanwise Variation of the Integrand of Equation (19) 
Displaying the Singularity at the Downwash Chord. 


In the following section, the singularity just described is identified as 
being logarithmic in its character and numerical evaluation of equation (19) 
will be accomplished using logarithmic quadrature integration functions. Also, 
the following section contains a description of how the spanwise integrand may 
be modified to provide a non-singular distribution that may be readily 
evaluated by Legendre-Gauss integration quadrature functions. 



IDENTIFICATION OF SPANWJSE SINGULARITIES 


The spanwise singularity contained within the integrand of the dipole 
downwash term was first identified by Multhopp (reference 8) in the develop- 
ment of the steady state lifting surface solution. The technique developed 
by Multhopp has been extended by Laschka (reference 9) to make it applicable 
for use in non-steady lifting surface solutions. 


The functional expression of the additional singularity contained within 
the dipole downwash term may be obtained by developing a Taylor series ex- 
pansion of AP(£,r))“ lkx 0 about the downwash station and performing the 
resulting chordwise integration. That is, let AP($,n)e~ lkx 0 be approximated 
by spanwise-chordwise Taylor series given by 


AP(S,n)e 0 = AP(x,y) + (n-y) |^[ AP (5.n) °] 

n=y 

S=x 



a r -i kxJI ~ 2 r -ikx" 


+ (5-x) 

|r AP(£.y)e ° + (n-y) a E JaP(«.n)e 0 



L - 5-x- “ n=yJ 


S=X 


The singularity in the spanwise integrand is identified by inserting 
equation (20) into the dipole chordwise integral of equation (10) and 
performing the indicated integration that results in the spanwise integral 
given as 



Thus, the additional singularity is identified as being logarithmic in 
character which is amenable to integration by use of logarithmic integration 
quadrature functions. 



The third integrals of equation (10) also contain a spanwise singularity 
that is proportional to -£nB 2 (y-n) 2 . The singularity is identified by insert- 
ing the Taylor series expansion of equation (20) into the third downwash term 
and performing the chordwise integration that results in a spanwise integral 
given as 




= ik I I AP(x,y)+(n-y) 


3 

9n 


AP(£,n)e °J 


£rt6 2 (y-n) : 


Regular 

Terms 


] dn 


( 22 ) 


n=y 


The fourth integral of equation (10) has a spanwise singularity that may 
be identified by redefining the logarithmic part of the kernel function to have 
the following form for 5<x 

-fnj((x-?) 2 +a 2 (y-n) 2 )^ - (x-£)J = £nB 2 (y- n) 2 - £n|((x-S) 2+ 3 2 (y-n) 2 )^ + (x-£)J 

Insertion of this expression into the chordwise integral results in the 
equation for I ^ of 

I 1 * 1 - - f - / [*>8’<y-n) 2 / AP(C,y)e' 1kX ° ♦ ^ularj d „ (23 ) 

_S X £ 

Once the singularities have been identified they may be subtracted out 
of the integrand and analytically evaluated outside of the integral. 

Figure 5 shows the very smooth spanwise distributions of the integrand 
I(y,n) obtained for the dipole term by subtracting out the identifiable 
singularities of equation (21). Removal of the singularities allows the use 
of Legendre-Gauss quadrature formulas to evaluate the spanwise integral using 
only a very small number of quadrature stations. As in any singularity sub- 
traction process, localized gradients may appear that may require additional 
subdivision of the integration interval to accurately evaluate the integral. 


It should be noted that there is a slope discontinuity yet existing at 
the downwash station and this may be smoothed by subtracting additional terms 
of the Taylor series expansion. 



Figure 5 — Spanwise Variation of Integrand of Equation (19) with Singularities 
of Equation (12) Removed from the Integrand 

As a consequence of the above there are several options open to the 
analyst in evaluating the spanwise portions of the downwash integral and 
are: 

(1) ignore the logarithmic singularity since it is of relatively low 
order > 

(2) subtract the singularities out of the expression and analytically 
evaluate their downwash contribution, and 

(3) apply numerical integration quadrature functions that can evaluate 
the logarithmic term at the downwash chord. 



Numerical investigations have been performed to determine the validity 
in using any one of the three integration options within the control surface 
program. The analysis planform shown in figure 6 is a rectangular wing having 
a full span flap deflected in steady flow. The deflection shape and step 
function downwash distribution are also displayed to show the discontinuity 
that must be matched by the mathematical downwash calculation of the integral 
equation. The assumed pressure distribution used in the analysis is composed 
of a chordwise singularity term x ) 2 and other modifying terms required 

to satisfy the boundary conditions along the planform edges. 



Figure 6 — Analysis Planform Used to Evaluate Effect of 

Spanwise Singularities on Downwash Distribution 

Results of applying option (1) to evaluate the spanwise integral of 
equation (19) are shown in figure 7. The integral is evaluated using the 
interdigiting process suggested in reference 7 where a single Gauss -Mehler 
quadrature function is applied over the complete span of the planform. The 
singularity at the downwash chord is not recognized by this integration 
procedure and consequently the resulting chordwise downwash distribution is 
smooth and continuous across the hingeline. It appears that lifting surface 
solutions obtained using option (1) will not provide a satisfactory mathe- 
matical representation of the boundary conditions for physical wing-control 
surface combination. 
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Figure 7 — Chordwise Downwash Distribution at Mid-Span Station Resulting from 
Applying Option (1) (logarithmic singularity ignored) 


Figure 8 presents the results using the second option where the 
singularity is treated by approximating the pressure terms by a Taylor 
series expansion and removing the singularity from the spanwise integrand. 
Equation (19) has been modified by subtracting the singularities (identi- 
fied in equation (21) ) and evaluating their effect outside of the integral. 



Figure 8 - Chordwise Downwash Distribution Resulting from Applying Option (2) 
(logarithmic singularity is subtracted and evaluated separately) 



The downwash distribution appears to be reasonable at large distances 
from the hingeline; however, the downwash becomes singular at the hingeline. 
This is to be expected since the pressure term used in the analysis to provide 
the step function change in boundary conditions across the hingeline is pro- 
portional to log |?-x c |, which cannot be approximated by a Taylor series in the 
vicinity of the singularity at 5 =x . 

V 

The downwash distribution shown in figure 8 indicates that the downwash 
collocation stations are to be restricted to chordwise locations that are 
spaced well away from hingeline where large downwash gradients are developed. 
This restriction may not be satisfied for analyses of wings having small per- 
cent chord control surfaces. 

This restriction may be removed by applying option (3) where the spanwise 
integral is evaluated by subdividing the integration interval, -s<n<s, into 
several subintervals and applying quadrature formulas appropriate to the 
functional characteristics of the integrand. 

Results of applying option (3) to evaluate equation (19) are shown in 
figure 9. The integration interval has been divided into subregions and 
appropriate quadrature formulas have been applied in the subregions, i.e., 
square root quadrature formulas applied at ends of the interval, logarithmic 
quadrature formulas used on either side of the downwash station, and Legendre 
quadratures applied within the remaining integration intervals. 



Figure 9. Chordwise Downwash Distribution Resulting from Applying Option (3) 
(logarithmic singularity evaluated by appropriate quadrature 
formulas) 



The resulting downwash distribution (figure 9) does contain the step 
function change across the hingeline necessary to satisfy the boundary condi- 
tions and does not exhibit any large gradients near the hingeline that may 
cause solution sensitivities for analysis of small chord control surfaces. 

Additional numerical investigations of calculated downwash distributions 
have been accomplished to identify the specific integration technique that 
would provide the least error in calculated downwash values for stations near 
the leading edge where large pressure gradients appear in the solution process. 
An assumed pressure function that is inversely proportional to the square root 
of the distance from the leading edge was applied to a rectangular planform 
having an aspect ratio of six. Figure 10 represents the steady state chord- 
wise downwash distributions obtained at the mid semispan station using option 
2 and option 3. The distributions appear to be similar for stations located 
at large distances from the leading edge, however, the option 2 distribution 



Figure 10 — Chordwise Downwash Distribution for Inverse Square Root Pressure 
Function Resulting from Applying Option 2 and Option 3. 



(where the logarithmic singularity is subtracted from the integrand and 
evaluated separately) is tending to become singular near the .leading edge. 
Again, this tendency for the calculated downwashes to become singular is 
expected since the Taylor series expansion of the inverse square root does 
not exist at the leading edge. The distribution obtained by applying 
appropriate quadrature formulas to the spanwise integrand (option 3) appears 
to be continuous for all regions including stations near the leading edge. 

Consequently, it appears that application of option 3 will provide 
solutions that will not be sensitive to the relative location of downwash 
collocation stations for analysis of either small percent chord control 
surface configurations or for analysis of standard lifting surface 
configurations. 

The above numerical investigations have been accomplished for the 
steady flow condition by evaluating equation (19) having the reduced 
frequency, k, set equal to zero. The following section establishes the 
numerical methods used to evaluate the downwashes for sinusoidal unsteady 
motion. 


Evaluation of the Unsteady Downwash Integral 

The downwash integral for unsteady motion is evaluated using a 
modification of the procedure suggested by Hsu (reference 7) where two 
singularities are subtracted and evaluated outside of the integral. The 
procedure previously described to evaluate the downwash due to the dipole 
term is modified by changing the definition of H(x,y,n) of equation (17) 
to become 


H(x,y,n) = • G(x,y ,y) j 

such that 


H(x,y,y) 


lim H(x,y,n) 
n-*y 


= a.G(.x , y,n) 

3n 


n=y 


G'(x,y,y) 


(24) 


(25) 
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Then equation (19) may be written in the revised form of 


B'^x.y.o) = -/(sV ) 4 6 (X,y 'fcyj X,y,y) ~ 6 '(»,y.o> 


n=y 


dn 

Tn=yT 


-s 


+ tt G(x,y,y)-G'(x,y,n) 


n’y /( s! V) % 


dn 

(n-y) 


( 26 ) 


-s 


and this is returned to another form using the definitions of AP(£,n) and 
the dipole kernel that results in 
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-ikx 
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ikx. 
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where AF(£>n) = AP(£,n)/(s 2 -n 2 )'* 
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The downwashes due to the non-singular part of the kernel function are 
evaluated using Gauss-Legendre integration quadrature functions. The final 
form used to evaluate the downwash integral is given as follows: 


s x t 

f (s 2 -n 2 ) % f APU,n)K(x,s»y»n) d(b 0 S)d(b 0 n) 


-s 


/(s 2 -n 2 ) % p(n)| 


C $ (y,n)+C ns (y,n)r+f(y) 


1 Zl* v v) s ' (x,y, ’ l) |n=yl 

.n)J + f(y ) ^ J 


dn 


+ tt G (x,y,y) + 7T y G'(x,y,n) 


n=y 


(28) 


where 


i p (5,n) = - f(n) g(5,n) 

(s 2 -n 2 )^ 

f(n) is the spanwise assumed pressure distribution divided 
by (s 2 -n 2 )^. 

g(?.n) is the assumed chordwise pressure distribution. 


C s (y,n) represents the chordwise integral due to the singular part 
of the kernel function 


C s (y,n) - 
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C n$ (y,rO represents the chordwise integral due to the non-singular 
part of the kernel function 


c„ s (y.o) - 



§(?»n) K(x,£»y>n)flg 


d? 


( 30 ) 


A 

G(x,y,y) = f(y) 2 J g(£,y) e' lk ( x "^ d£ 


(31) 


( x t 

f(n) / g(5,n)e~ lk(x ~ c) [ 1 + . X ~S -1 dg 1 (32) 

{ L V( x ^) 2+ e 2 (y-n) 2 J ] 

1 n=y 


The chordwise integrals of equations (29), (30), and (31) are numerically 
evaluated by subdividing the integration interval into small regions and apply- 
ing integration quadrature formulas appropriate to the functional character- 
istics of the integrand. The spanwise integration interval of equation (28) 
is also subdivided into small integration intervals that includes an interval 
around the downwash chord where logarithmic integration quadrature formulas 
are applied to evaluate the spanwise singularity previously identified. 

Evaluation of the spanwise derivative defined by equation (32) may be 
accomplished following the rule for differentiating a definite integral 
providing the integrand g(£,n) does not contain terms that are singular. 
Special procedures are developed to evaluate equation (32) for the case of 
having pressure functions that have an inverse square root singularity at 
the leading edge or having a logarithmic singularity at the hingeline and 
are given as follows. 
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The chordwise Integral for the case of having a pressure function that 
has an inverse square root singularity at a swept leading edge is given by 


I = f(n) f * -===r gU,n)e“ ik(x "^ Fl + ^ = =1 

' yj s-x^ [_ )/(x-£) +B (y-n ) 2 J 




(33) 
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It is required to evaluate 


n=y 


e) 2 +e 2 (y- 

and may be accomplished by forming a 


Taylor series expansion of the last term of the integrand about n=y; that 
is 
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and apply a transformation to reduce the integrand to an easily integrated 
form by 1 etti ng 


T “ ? 


then 


d? = dx 


V x t 


1 • f(n)e ik *fc / g(T + x r n)e- fk(x ' T > 1 
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e- ikXk f X ^ 9 (^,„,e- 
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The last integral of the above expression will not contribute to the 
derivative value when evaluated at n=y and may be neglected in the 
derivative calculation. 



Following the rule for differentiating a definite integral and applying 
the chain rule 

9F(U(n) »ri) _3f.9U.3F 
9n " 9U 9n 9n 


and noting that 

3 1 1 *m] 
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where 

6 (z) is the Dirac delta function. 

then the expression for the spanwise derivative becomes upon applying the 
inverse transform 
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A similar expression may be obtained for a chordwlse Integral that 
contains a logarithmic singularity such as 


I = f(n) J £o3(€-x ) 2 g(€,n)e~ ik t*~ g) Tl + — .1 

i L V( x -o*+* a ( y-n) a J 


d? (37) 


where x £ is a hingeline that lies within the integration Interval. 


The integration interval is divided into two parts 


t 
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lim 

e-K) 


x c" e x t 
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x c+e 


and the same procedures are applied resulting in the expression for the 
spanwise derivative given as 


31 

an 


= Mini 


n=y 


3n 


n=y 


A 

2 f £<?gU-x c ) 2 gU,y)e’ lk ( x_ ^ d£ 


+ f(y) 2 J log ( g-x c ) 2 + MiiUl j e -ik(x-5) ^ 
y n y 

X 

+ f (y) 2ik J log(Z-x c ) 2 g(£;,y)e~ lk (*~Z) d? 


3 x 

- f (y) 2 3 ^£og(x-x c ) 2 g(x,y) 


(38) 


31 



Thus, the spanwise derivative of equation (32) may be evaluated for func- 
tions that have singular or non-singular characteristics and the downwash in- 
tegral of equation (28) may be evaluated using integration quadrature functions 
that are appropriate to the functional characteristics of the integrand. 

Lifting surface downwash distributions may be generated by evaluating equa- 
tion (28) for any pressure expression that may be either continuous or non- 
continuous over the surface. Solutions for configurations having discontinuous 
downwash distributions require application of special techniques to minimize 
the solution sensitivity. 

Solution Process for Discontinuous Downwash Distributions 

No direct means is available to obtain a closed form solution of the down- 
wash integral equation that satisfies the boundary conditions exactly everywhere 
on the surface. Approximate solutions are used in lieu of an exact solution 
process. The approximate solution process is one of generating a finite set of 
downwash sheets and constructing linear combinations of these sheets to satisfy 
the boundary conditions of a finite set of pre-selected control points on the 
surface. 

The downwash sheets are obtained through an evaluation of the downwash 
integral using assumed pressure distributions that satisfy the required loading 
conditions on the planform edges. The assumption usually made is that if the 
boundary conditions are satisfied at a suitable number of control points dis- 
tributed over the surface, then boundary conditions over the rest of the sur- 
face will be satisfied within small error limits. 

This approximate procedure is equivalent to assuming that the downwash 
sheets may be represented by a series of polynomials that are continuous over 
the surface and may be combined to satisfy the surface boundary conditions at 
arbitrary control points. 
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This procedure provides reasonable results for analyses of configurations 
having continuous downwash distribution; however, it may fail to produce con- 
sistently accurate results when applied to analyses of wing-control surface 
configurations that have discontinuous downwash distributions. 


The solution procedure used for analysis of continuous downwash configura- 
tions has been modified such that solutions may be obtained for configurations 
having discontinuous downwash distributions and is described in the matrix 
equation 


|W(x,y) 


-| // 


APU,n) ae K U*n) d5dn 


•where 


[// 


AP-U.n) KU.n) d£dn 


x K 


|W(x,y)| is the kinematic downwash obtained from the motion 

of the lifting surface and containing discontinui- 
ties along the edges of a control surface. 


J*J*APU.n) ae KU>n) dedn 


is the surface downwash distribution that is 
formulated having identical discontinuities along 
the control surface edges as defined for the 
kinematic distribution. 


AP(5»n) a _ is the pressure distribution developed by the 

Qc 

asymptotic expansion process that will provide the 
required change in boundary conditions across the 
control surface edges defined by the kinematic 
distribution. 


yyy.U,n) KU,n)d£dn 


is the standard lifting surface downwash distribu- 
tion that is smooth and continuous and is calcu- 
lated using smooth and continuous assumed pressure 
distributions AP.(g,n) 
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The solution process is a two step process. The first step is to calcu- 
late the "residual downwash distribution by performing the operation 


W(x,y) 

kinmcuUc 



K(e,n) d^dn 



W(x,y) 

KulduaJt 


The residual distribution is to be smooth and continuous and is used as 
the boundary condition for the second step, which is the solution of 


W(x,y) 

neAlduaZ 




U.n) K(c,n) dedn x a A 


The final pressure distribution is then the sum of the pressure compon- 
ents, i.e. 

APU.n) * p (£.n) ae + £ AP^(C,n) 

i 

A schematic of the solution process is presented in figure 11 for a con- 
figuration having an oscillating trailing edge control surface. 
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The kinematic downwash distribution is developed from the surface motion defi- 
nition and displays discontinuities along the boundaries of the control surface. 
A downwash distribution is calculated having identical discontinuity values 
that are contained in the kinematic distribution and is obtained by evaluating 
the downwash integral of equation (28) using an asymptotic pressure expression 
developed for a swept hingeline trailing edge control surface. The residual 
distribution is then formed by subtracting the kinematic distribution from the 
calculated distribution. The residual distribution must be smooth and con- 
tinuous and is used as the new boundary condition definition for which solu- 
tions are obtained using standard lifting surface procedures. 

A similar procedure is used to obtain solutions for motions of a small 
chord leading edge control surface as is shown in figure 12. 




QOmiMSH - FULL CHOW LEAPING 
EPGE CONTROL SURFACE 


POWNUASH - LARGE CHOW TRAILING 
EPGE CONTROL SURFACE 



RESIPUAL POVHUASH 


Figure 12. Solution Process for Leading Edge Control Surface 
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The residual distribution is formed by adding to the kinematic distribu- 
tion a downwash distribution of a large chord trailing edge control surface 
that has a hingeline located at the same position as the small chord leading 
edge control surface and then subtracting out the downwash distribution gener- 
ated for a full chord control surface. 

The key to this solution process is the development of pressure distribu- 
tions having coefficients of known strengths that will provide identical values 
of downwash discontinuities that are contained in the kinematic downwash dis- 
tribution. Formulation of the analytical pressure distributions that provide 
the required discontinuities in boundary conditions are briefly discussed in 
the following section. 


Pressure Functions 

Use of the modified solution process requires two sets of pressure distri- 
butions to be defined over the lifting surface planform. One of the distribu- 
tions is to provide identical values of downwash discontinuities as defined by 
the kinematic downwash distribution and the second pressure function is re- 
quired to satisfy the smooth and continuous boundary conditions defined by the 
residual downwash description. 

Incremental pressure expressions have been developed by Landahl 3 to satis- 
fy the change in boundary conditions around the edges of a non-swept hingeline 
trailing edge control surface configuration that has its side edges inboard of 
the planform tip. 

The solution of Landahl 3 was obtained using the asymptotic expansion pro- 
cess. The behavior of the pressure functions is readily obtained for regions 
in the immediate vicinity of the control surface leading edge and side edge 
since the effect of using asymptotic expansions is to magnify the singular 
region by stretching the coordinates by a small quantity e . 
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The solution derived by Landahl 3 that will exactly satisfy the required 
change in boundary conditions near the corner of a non-swept hingeline 
trailing edge control surface is given as 


m.two) * - 

. Mill 


[l+2ik(£-£ c )] log | ^(S-E c ) 2 +B 2 (n-y s ) 2 - 6(n-y $ ) 


( 39 ) 


^nl[2ik-k 2 (?-c c )] (n-y s kog|y(?-e c ) 2 +6 2 (n-y s ) 2 - (5-€ c ) 


where 0(n) is the rotation angle of the control surface measured in a plane 
perpendicular to the y axis. 


The expression contains the expected logarithmic singularity along the 
hingeline and the spanwise loading has slope singularity at the side edge of 
the control surface that is known to exist for lifting surface loadings on 
configurations having discontinuous downwash distributions. 

Pressure functions have been formulated by Ashley 10 for a swept hinge- 
line trailing edge control surface that extends to the tip of the planform. 
Application of asymptotic expansions along with use of Fourier transforms 
results in the following first order solution given in reference 10 and 
expressed in integral form as 



where the bar coordinates represent stretched coordinates measured from 
hingeline tip station, l.z., x = (?-x c )/e, y = (n-y $ )/e 


The integrals may be expressed in a more easily handled integrable 
form by making the substitution 
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and performing the indicated integrations which result in the shortened 
expression 
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where 
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The above expression appears to be correct since Cp^(x,y,0) 
approaches zero having a slope that is proportional to the square root of the 
distance from the wing- tip and also contains a term that is proportional to 
£og(s;-x c ) 2 to provide the proper streamwise change in boundary conditions 
across the hingeline. 

Pressure expressions developed in appendix 'A' satisfy the boundary 
conditions around the edges of a arbitrary configuration having a swept 
hingeline trailing edge control surface with side edges inboard of the 
planform tip. Asymptotic expansion procedures are used to reduce the very 
complicated boundary value problem to three simpler boundary value problems 
for which solutions are obtained by making a direct application of Green's 
Theorem. The resulting pressure expression obtained from the solution of the 
three boundary value problems is given as ( with x * $-x and y = n-y* ) 

V J • 
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The expression for F(x,y,0) may be simplified by expanding the 
exponential, carrying out the indicated multiplication and collecting terms 
that are linear in x and y which results in the expression 


P(x,y,0) 


| [l + 2ik(5-C c ) + i^<e-5 c )] 

J-pn - k ! (s-? c )J y loq f V x 2 +B ! y ! - xj (43) 


The limiting value of this expression for A=0 is equivalent to the 
expression developed by Landahl 3 for the rectangular non-swept control 
surface. The above expression for A=0 is 


P(x,y,0) 


i [l + 2ik(£-e c ) + tog [V(^ c ) 2 +B 2 y 2 - Sy] 

| f2ik-k 2 (5-e c )J y ^og^V(C-C c ) 2 +B 2 y 2 - (€-e c )J ( 44 ) 


and is identical to Landahl's zero sweep expression. The expression developed 
for the arbitrary sweep angle case does exactly satisfy the change in bound- 
ary conditions across the hingeline and side edge of a swept hingeline trail- 


ing edge control surface and should provide smooth and continuous residual 
downwashes for which solutions may be obtained that are relatively insensitive 
to the number of collocation stations placed on the lifting surface. 

Pressure expressions have also been developed in Appendix "B" that satis- 
fy the change in boundary conditions across the side edge of a full chord 
swept leading edge control surface. 

The zero-th order solution is given in terms of an inverse square root 
function that is multiplied by a coefficient whose value can only be deter- 
mined after performing a global integration. The inverse square root function 
is continuous in the spanwise direction, that is, it does not exhibit any 
spanwise singularities that could contribute to the change in downwash across 
the side edge. Consequently, the zero-th order boundary value problem solu- 
tion is not included within the incremental pressure expression necessary to 
satisfy the change in boundary conditions across the control surface side edge. 

The pressure expressions obtained as solutions of the first order and 
second order boundary value problems are given as follows: 

First order solution is: 
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The second order solution is: 
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where Cg ■ v ; v ■ e H (ik$) ; x * e H (l-ikx c ) 

3 

X - 5-x ls ; y • n-y s 

( note: x,y used here is not to be confused with 
the downwash coordinate x,y ) 

Pressure distributions obtained by the above described asymptotic ex- 
pansion process are valid only in localized regions of the planform since the 
distributions of the outer regions have not been defined nor matched with 
the inner region solutions. However, the essential part of the pressure 
distribution that provides the discontinuities in downwash along the control 


surface edge is defined by the inner region solution. The distribution in 
the outer regions may take on any convenient form that satisfies the required 
loading condition along the edges of the planform since a unique distribution 
is not being determined at this stage of the solution process. Consequently, 
the surface pressure loading functions that provide the downwash discontinu- 
ities are obtained by extending the applicable range of the inner solution 
and modifying these distributions to meet the condition of having the pressure 
vanish in proportion to the square root of the distance from the edges. 

The functions used to modify the asymptotic pressure expressions are 
usually artfully devised and are not unique expressions. However, there are 
certain characteristics that are required to be satisfied by these functions 
such that the resulting residual downwash sheet will contain a minimum amount 
of distortion. Some of the required functional characteristics are established 
in the following section. 

Pressure Modification Functions 

Pressure loading functions are developed using the asymptotic pressure 
expressions and extending their applicable range from localized regions to 
cover the entire planform. 

For the case of a trailing edge control surface, the asymptotic pressure 
expressions that are valid only for small regions around the discontinuity 
region may be written as 
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_ Jh_ C 2 log / U-x cs ) 2 +e 2 (n-y s ) 2 - (e-x cs ) 


( 47 ) 


where 


C x « CjA-ej/wB) 


Cj , Cg are defined in equation (42) 



C(x,y,0) Is related to the change In pressure loading across the 


surface by 


ac p - * 2 S, 


AP _ . „ „ 

- + 2«p t 


AP = + C n x pV 2 
HI 

where C is given by the above equation (47) 


Then 

&P(e*n) - pV 2 


t eu _ r / r — - 3 2 (n-y s )+(c-x JtanA "I 

~C X tog U-x cs ) 2 +e 2 ( n-y s ) 2 £ J 


®H — 

— C 9 log / (C- 
ir ^ L 


x cs ) 2 +6Mn-y s ) 2 - U-x cs ) 


]) 
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A characteristic of the first logarithmic term is that it is singular 
only along the hingeline and the coefficient value of the singularity must 
be maintained along the hingeline such that the unit change in downwash is 
obtained for all points on the control surface. Thus, to distribute this 
function over the entire lifting surface a pressure modifying function is 
devised to multiply this term such that the coefficient value is maintained 
along the hingeline and approaches zero at the planform edges in proportion 
to the square root of the distance from the planform edges. 

The modifying function is made up as a product of two functions such 
that the streamwise and spanwise boundary conditions may be satisfied 
independently. 


The spanwlse modifying function, H(n), is defined as 


{ i \fnr (1. + . 5C + .375C 2 ) In! > y 0 

(49) 

1.0 Ini < y 0 

g _ 

where L (s-y Q ) 

y Q = the spanwise coordinate of outboard side edge of the 
control surface. 

H(n) has a value of 1.0 and a slope of zero at the control surface outboard 
side edge and then diminishes to zero in proportion to the square root of the 
distance from the planform tip. 

The chordwise modification function, Ej, must be devised such that Ej = 

1.0 at the chordwise hingeline station to maintain the singularity strength 
and must approach zero in proportion to the square root of the distance from 
the leading and trailing edges. 

However, there is an additional condition required of Ej such that smooth 
and continuous residual downwash are developed. The additional condition im- 
posed on Ej is that E^ must have a zero slope in the streamwise direction at 
the hingeline station. The requirement of having a zero slope may be illus- 
trated by examining the results of a numerical investigation that was accom- 
plished to identify the essential properties of Ej such that reasonable distri- 
butions of residual downwashes would result. 

Downwash distributions were obtained at the four downwash chords shown in 
figure 13 by applying the zero sweep loading function given as follows: 

AP(e,n) * 4pV 2 f(n) H(n ) tog j/ U-e c )2+62(n-y s ) 2 - e(n-y $ )] (50) 



Y/S - SVAM1SE COORDINATE 

Figure 13. Analysis Plan-form Used to Evaluate Ej Characteristics 

One of the many modifying functions investigated is the function suggested 
by Landahl 3 which has a value of 1.0 at the hingeline and approaches zero in 
proportion to the square root of the distance from the leading and trailing 
edges as shown in figure 14. 



o .2 .4 .6 .1 1.0 


x/c- CHomnsE coordinate 

Figure 14. Functional Characteristics of E^U) 
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The residual downwashes that result using the EjU) function of figure 14 
are shown in figure 15 and represent the distributions on the aft portion of 
the wing-control surface regions. The residual downwashes are smooth and con- 
tinuous for downwash chords that are inboard of the control surface side edge. 
However, for spanwise downwash stations located on the control surface, the 
distributions contain a slope discontinuity that may not be easily satisfied 
by the downwashes generated by the remaining portion of the solution process. 



Figure 15. Residual Downwash Distributions Resulting from Using the 
EjU) Function of Figure 14. 

This slope discontinuity could have been anticipated by noting that the 
logarithmic term of the loading function may be written as 

log(t-Z c ) 2 - 4og[/ U-e c ) 2+ B 2 (n-y s ) 2 + B(n-y s )] (51) 
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and coupling this with the E^(c) function represented as a Taylor series ex 
pansion about the hingeline in the form of 


E l"> ■ 1 ♦ «-'e> f>| + (i r c>2 0 1 


(52) 




then insertion of these terms into downwash equations provides a series of 
downwash integrals having the properties presented in figure 16. 


JJfU.n) £ogU-5 c ) 2 K(x,5.y»n) d?dn 



fU.n) (e-C c ) £ogU-e c ) 2 K(x.5,y,n) d5dn 


JJfU.n) U-e c ) 2 £ogU-e c ) 2 K(x,c,y, n) d^dn 



Figure 16. Chordwise Downwash Variations Caused by 
Singularity Multiplier Functions 

It is apparent that the finite angle slope discontinuity displayed in 
figure 15 is due to the non-zero slope at the hingeline contained in the 
definition of E^s) shown in figure 14. Also, it may be noted that large 
curvatures in downwashes may be minimized by devising a function that con- 
tains only a small amount of curvature in the vicinity of the hingeline. 

After extensive numerical evaluations of various EjU) functions, the 
function shown in figure 17 represents the final form of E^(?) that has been 
selected to modify the chordwise loading expressions. 




Figure 17- Final Form of Chordwise Modification Function 

Residual downwash distributions of figure 18 were obtained using the 
EjU) function of figure 17. It appears that this function does provide slope 
continuity across the hingeline as well as developing minimum waviness over the 
rest of the chord. Solutions using this modification function should then be 
relatively insensitive to the distribution of downwash collocation stations. 



Figure 18- Residual Downwashes Obtained Using EjU) of Figure 17 
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The second logarithmic term of equation (42) contains a singularity of 
the type (n-y ) log (n-y ) that is to be applied along the side edge of the 
control surface. The strength of this singularity is to be maintained along 
the length of the side edge to provide the proper downwash discontinuity along 
the side edge for the unsteady flow condition. 

A function (suggested by Landahl 3 ) that may be used as a pressure modifi- 
cation function to maintain the singularity strength along the side edge as 
well as forcing the pressures to zero in proportion to the square root of the 
distance from the trailing edge is given by the following equation and also 
presented in figure 19. 



Figure 19. Schematic of n)» Pressure Modification Function 
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The side edge modification function given by equation (53) has been 
investigated using numerical procedures to determine its range of validity in 
maintaining the side edge singularity strengths. Results of the investigation 
indicate that the function requires a modification to reduce its spanwise 
gradients in the region of the side edge-trailing edge junction. It appears 
that in applying equation (53) to the side edge pressures that the singularity 
strengths t ._ar the trailing edge are being modified to such an extent that the 
spanwise integration quadrature functions are unable to detect the infinite 
slope of the pressure function which results in an error in predicting the 
spanwise downwash discontinuity. 


Equation (53) has been modified to provide a proper evaluation of the 
side edge downwash discontinuity and is given in the following form for use in 
a full chord control surface analysis. 


^2^»n) 


<v 5)2 

( x t - x *) 2+ C(n-y s ) 2 


<V\> 2 

2 - 2 
(x t -e) + c(n-y $ ) 



(54) 


Equation (54) has a value of 1.0 along the side edge and leading edge and 
also forces the pressures to approach zero in proportion to the square root of 
the distance from the trailing edge. 


The side edge modification function used for a trailing edge control sur- 
face configuration is given as 




U-x t ) 2 

% 


(c-x Jl ) 2 +C(n-y s ) 

• 

(x t -0 + C(n-y s ) 


(55) 


In this case, the pressures are forced to zero in proportion to the squan 
root of the distance from both the leading and trailing edges. 
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Loading Functions for an Inboard Partial 
Span Control Surface 


Pressure loading functions formulated for a trailing edge partial span 
control surface having side edges Inboard of the planform tip are devised 
from the trailing edge control surface pressure expansion of equation (42) as 
follows: 



e H is the control surface rotation angle relative to the wing 
rotation angle and measured in a plane perpendicular to the 
y axis. 

SF is the symmetrical sign factor that has a value of SF = 1.0 
for symmetrical loadings and SF = - 1.0 for anti symmetrical 
loadings. 

D 

The g^U.n) loading functions are the loading functions that are due to 
the right-hand semi span control surface deflection. The functions are formu- 
lated such that the hingeline singularity is turned on at the inboard side 
edge and maintained at full strength along the hingeline and then turned off 
at the outboard side edge. 

D 

The gj(c, reloading functions are non-singular for spanwise stations 
that are outside of the spanwise control surface location. Consequently, 
the loadings are extended over the rest of the surface to the planform tips 
where the loadings are forced to zero in proportion to the square root of 
the distance from the tips. 
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gju.n) » Ej H(n) Cj ^to§ [b/(C-x c1 ) 2 + ^(n-y^ 2 - (6 2 (n-y i ) + U-* ci )tanA c ) ] 
- &*[? /(c-x c0 ) 2 + B 2 (n-y 0 ) 2 - U 2 (n-y 0 ) + (?-x co )tanA c )] J (57) 


The subscripts i and 0 represent the inboard and outboard side edges 
respectively. The parameters , y.. and x Q , y Q represent the coordinates of 
the hingeline at the inboard and outboard side edge stations of the control 
surface. 


The Ej function forces the loadings to go to zero in proportion to the 
square root of the distance from the leading and trailing edges and also 
has a value of Ej = 1.0 and zero streamwise slope at the hingeline. 


The expression for E^ for stations ahead of the hingeline is given by 


/U-x t )(2x c -x £ -g) 

( V x t ) 


( 58 ) 


for stations aft of the hingeline, 

/(x t -0(C-2x c +x t ) 

E i (vg — 


( 59 ) 


The H(n) function is the spanwise pressure modification function having 
a value of H(n) =1.0 for spanwise stations located between the left hand 
outboard and right hand outboard control surface side edges. For stations 
outboard of the outboard side edges, the H(n) function goes to zero in pro- 
portion to the distance from the planfonn tip station: 


I 


52 



/l. - l (1 + . 5C + .375C 2 ) 


(60) 


H(n) * 


1. 


V*here 


. Clnl - y 0 ) 

c » — 

(s - y 0 ) 


Ini iy 0 
|n| < y 0 


The Cj function represents the coefficient multipliers obtained in 
deriving the pressure expressions by means of asymptotic expansions. The 
Cj expression of equation (42) has been shortened by expanding the exponen- 
tial function and performing the indicated multiplication that results in 
the expression 


-L 1 + 2ik(5-C c ) +^pU-C c ) - U-n c ) 2 ( 
3 L 


k 2 M 4 A k 2 M 2 + 


4S 2 8 2 B 2 


-)1 

2 / ] 


(61) 


where $ ■ + tan 2 A c 


The loading functions that are due to motions of the left hand side 
control surface are obtained from the above expression by replacing y Q , y.. , 
and tanA by -y , -y. and -tanA , respectively. 

w U I v 

The second part of the pressure function defined in equation (42) is 
formulated into loading functions g^U.n) and g 2 U»n) by applying other 
pressure modification functions. The loading function due to the singu- 
larity developed for the right hand side control surface motion is given as 

gJU.n) - H(n) £ E 2 (y t ) C^) ( n - y i ) log [/(t-x c1 ) 2+ B a (n-y i ) 2 
- E 2 (y 0 ) C 2 (y 0 ) (n-y 0 ) log [/(c-x c0 ) 2 ^ B 2 (n-y 0 ) 2 
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The subscripts i and o pertain to the inboard and outboard side edges 
respectively. The (x^, y..) t (x Q , y ) parameters are the hingeline coordinates 
at the inboard and outboard side edges respectively. 

The logarithmic function is singular only along the side edge which 
requires that the strength of the singularity be maintained only along the 
side edge and then may take on any other desired form away from the side 
edge. The E 2 (y s ) ^ unct10n 1S devised in such a manner as to have a value 
of E 2 (y s ) = 1.0 along the side edge to maintain the singularity strength 
and is also devised to have a characteristic of reducing the pressures to 
zero in proportion to the square root of the distance from the leading and 
trailing edges. 


Using y $ to represent either y^ for inboard or y Q for outboard side 
edges results in the following chordwise pressure modification function. 


W 


(x t -0 2 

* 

r 

CM 

X 

1 

ur 

1 

_(x t -5) 2 + C( n -y $ ) 2 _ 


_U-X £ ) 2 + C(n-y s ) 2 J 


(63) 


where C * 0.10 


The spanwise modification function H(n) is the same function previously 
described and is given by equation (60). 

The C 2 (y s ) terms are obtained from equation (42) by expanding the 
exponential function and performing the indicated multiplication that 
results in the expression: 

c 2 (y$) * [ " 2ik + k2 (S‘* c ) + _L 2 (n-y $ )tanA c J (64) 

The loading functions developed for motions of the left hand side 
control surface are obtained from the right hand side terms by replacing 
y 2 » y Q > and tanA c by -y 2> -y Q , and -tanA c , respectively. 
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Loading Functions for a Full Chord Control Surface 


Loading functions for a full chord control surface are formulated in 
a similar manner using the pressure expressions of equation (45) and equa- 
tion (46). 


Loading functions developed from equation (45) for the inboard side 
edge of the right hand side control surface are given in localized coordin- 
ates of x » C-x t1 and y « n-y 1 . 


•Ju.,1 - e ,kH2;/B2 H(n.y,) 
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A - ytanA t 
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;[c 6 [v'x 2 + B 2 y 2 + x] + Cy [»6 2 + e 2 y 2 + x] J 


+ E 2 C 9 


_By_ yf t S «*(x - ytanA^I [/x 2 ♦ g 2 y 2 4 
|y| 2 A 2 + B 2 y 2 - x +ytanA t 


(65) 
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The spanwise modification function H(n.y.j) is defined as having a value 
of 1.0 with a zero slope in the spanwise direction at the control surface 
side edge and then diminishes to zero in proportion to the square root of the 
distance from either planform tip and is given by the following expression: 


H(n.y $ ) * /l - C [l + .5C + .375C 2 ] 


( 66 ) 


where 

C 


-(n - y s )/(|y s l + *) 

_ (n - y s )/( S - |y $ |) 


if n < y s 
if n > y s 


The modification function Eg has a value of 1.0 along the side edge and 
diminishes to zero in proportion to the square root of the distance from the 
trailing edge and is defined as 


E 


2 


(x t - C) 2 [(x t -x t ) 2 ♦ -Hn - y s ) 2 ] 1 * 
(x t - x 1 ) 2 [(x t - C) 2 + .l(n - y s ) 2 ] _ 


(67) 


The E^ modification function has a value of 1.0 along the leading edge 
and approaches zero in proportion to the square root of the distance from 
the trailing edge and is given by 

E 3 -[/ 3 - 2§ - § 2 ]/ 2 (68) 

where £ Is defined within the expression of 
{ * <«t + 5 t )/2 + [ (? t • E i )/2 ] 5 
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The pressure coefficients Cg through Cg are obtained from the expressions 
of equations (45) and (46), and are given as follows: 


c 3 - 2^r[i - iMx cs - * u )]/ i t 2 

C 4 ■ 2^[l - ik(x cs - x ts )](<k/8 2 )(l/», 2 ) 

C 5 - 1k[l - 1k(x cs - x ls )]/ B, 2 

O'. 24 


f faxtanA -(3tan 2 A -e 2 )yT 

ltanA t + (ik/$ 2 ) 1 - - J 


: 7 * 2^2 k 
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[l ♦ (1k/U 2 )(x - ytanA 4 )/2j 


-k 2 M 2 (x -x t ) + 1(B 2 - M 2 ) 


Cg « k' ^ 15 


j|l + (1k/e 2 )(x - .75ytanA £ )J 


C 9 • (1k 2 /2B 2 )£ 


-k 2 M 2 ( x c S - x ts ) + 1(e 2 - H 2 ) 
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The loading function described in equation (65) is formulated using pres- 
sure expressions that exactly satisfy the boundary conditions at the Inboard 
side edge of the full chord control surface located within the right hand semi 
span region. Loading functions due to downwash discontinuities at other side 
edges are obtained from equation (65) by inserting new definitions of x and y 
into the equation such that the equation is transformed into local coordinates 
of Individual side edge stations. 


Loading Functions Used in Solution of Residual Downwash Distribution 


The loading functions used to obtain solutions of the residual down- 
wash distribution are defined as follows: 

APU.n) • 4pV 2 /s 2 - n 2 £ 2 a n m g* m) U.n) (69) 

s n*l m»l 


where « nm are the unknown coefficient multipliers 


f (n) (n) * 


sin(n 4 >) 

sin* 


[ 1,3,5,-- 2n-l — 2N-1 symmetrical 

_ 

2,4,6,— 2n 2N anti symmetrical 

n » 1,2,3, -N 

N is the number of downwash chords on a semi span 


* 


cos _1 (n/s) 


9 (m) U,n) 


"cot ^ 

_ sin(m-l)$ 


m * 1 

m « 2,3,4, — M 


M is the number of downwash stations on a downwash chord 


The chordwise coordinate e is not to be confused with the control surface 
rotation angle © H ( n ) • 
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The loading functions g(s,n) have been described only in terms of a 
particular side edge deflection to indicate how the asymptotic pressures have 



been modified and extended to the planform boundaries. The complete set of 
expressions are provided in the computer program and are well annotated for 
each of the control surface configurations and need not be repeated in this 
section. 

Output of the solution portion of the program using the above described 
loading functions results in generalized forces described in the following 
section. 


Generalized Forces 

Generalized forces are formulated using two sets of pressure distributions 
that result from the solution of lifting surfaces having discontinuous down- 
wash distributions. 

Generalized forces are defined as follows: 

Q ij ‘ // AP j (c * n) H i U * n) d (V> d(b 0 n) (70) 

&uA&a ce 

where the limits of integration are the physical boundaries of the half span 
configuration. 

H.j U,n) is i th deflection mode shape. 

AP.(s»o) is the sum of the singular (asymptotic pressures) and regular 

J 

pressures obtained as a solution of the residual downwash 
distributions . 

APjU,n) - APj(c,n) r + APj(s,n) ae (71) 

The dimensions of AP^U.n) are: (force/area) (per unit q^) 
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The tegular pressure distribution APj(5»n) r is given as 


/s 2 - n 2 


N N 


4Pj{C. n ) r * 4,V 2 £ Z ‘m f M M 9 M U.«) 


( 72 ) 


n«l m-1 

w *' ere a nm^ ®re the coefficient multipliers obtained from the 

Integral equation solution using the residual downwash 
distributions as boundary conditions for the jth mode. 


f^(n) represents the assumed form of the spanwlse pressure 
variation defined in equation (69). 

g< m >U.n) Is the assumed form of the chordwise pressure 
variation defined in equation (69). 


The singular asymptotic pressure distribution Is given as 

4 

4P < (e,n), fl * 4pV 2 /s 2 - n 2 21 A. i f(i) 9(6»n) ( 73 ^ 

3 86 k-1 ' 

where f(n) * l/[4**4 2 - n 2 ] 

g(6,n) » are the loading functions that are formulated using the 
asymptotic pressure expressions and the pressure modifi- 
cation functions given in equation (56) or equation (65) 
depending upon the configuration being evaluated. 

(j) 

\-l are f° ur coefficients of a cubic polynomial used to 
define the spanwise twist distribution of e^(n). 



The program output of pressure is lifting pressure coefficient per unit 
generalized coordinate q. , and its units are: per unit generalized co- 

J 

ordinate q^ . Its numerical values are proportional to the modal -displacement 
input. The program output of pressures is given as follows: 


aP j ( K »n) 

* PROGRAM OUTPUT 

% pV 2 


( 74 ) 


where 


PROGRAM OUTPUT * 



N M 

Z Z •£> f (n) (n) * 

n*l m«l 


4 (j) 

2 \-i g(C.n) 
k-1 


Sectional generalized forces may be obtained for any number of arbi- 
trary locations along the semi span and are defined as: 

b 0*t 

q?j * hp»* j [ i ir?] H ' u - n)d(b '> ?) c»> 

Vt 


The section generalized forces Q^. have the units: (force per unit 

planform dimension per unit q.)*(units of modal -displacement input). 

J 

The program output of section generalized forces has the units: (plan- 

form length uni ts)x(modal -displacement units) per unit q. , and the program 
output is proportional to the modal displacement squared and to planform 
length. 
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The program output of section generalized forces is given by 

s 

Qjj - PROGRAM OUTPUT (76) 

>s pV 2 

where PROGRAM OUTPUT - 8 i& 2 - n 2 v 

s * 


Vt ~ ~ 

J[z Z. 

Ln*l m*l 

Vt 


nm 


f (n) M 


,(m) 


< 1 

U»n) + £ A k_i f (i) g(C.n) H i (^,n) d(b Q 0 

k*l ' J 


Total generalized forces are formulated using the definition of equa- 
tion (70) and the generalized force elements Q. . have the units: 

* J 

(force per unit q. )x(units of modal -displacement input) 

The program output of total generalized forces is given by 

q V Vt 

— !L • PROGRAM OUTPUT * 8 f f H.U.n) - n* x 

■s p V 2 J J 1 s 

0 Vt 

[ N M 4 (j) "| 

2 X + 2 A k-1 9(C»n) d(b Q c) d(b Q n) 

n»l m*l k*l ~ -» 

(77) 

The program output has the units: (planform area units)x(modal dis- 

placement units) per unit q. , and the program output is proportional to 

J 

the modal displacement squared and to planform area. 
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It should be noted that the sectional generalized forces of equation (75) 
are used to obtain information on the spanwise loadings. There is no restric- 
tion on the number of stations where sectional generalized forces may be ob- 
tained and is limited only by user option. 

The spanwise locations where the sectional generalized forces are calcu- 
lated do not coincide with the spanwise quadrature stations where chordwise 
integrals of equation (77) are evaluated. The chordwise integrals of equa- 
tion (77) are calculated at particular spanwise stations depending upon the 
quadrature distribution being used for a particular configuration and are not 
printed as output from the program. 

PROGRAM CAPABILITIES AND LIMITATIONS 

The computer program developed around the preferred solution process may 
be used for a vast array of analysis configurations and has capabilities and 
limitations listed below: 

1. Trailing edge control surface configurations may be composed of a full 
span, inboard partial -span, outboard tip partial -span, or up to six 
individual control surfaces including side-by-side common edge control 
surfaces, and nested controls. 

2. For reasons of detail accuracy, leading edge control surfaces are recom- 
mended to be limited to control surfaces having side edges inboard of 
the planform tip station. 

3. All hinge lines of trailing edge control surfaces are assumed to be 
located at the leading edge of the control surfaces. All hinge lines of 
leading edge control surfaces are assumed to be located at the trailing 
edge of the control surface. Control surface hinge lines may have indi- 
vidual definitions; however, each hinge line is assumed to be defined by 
a linear equation. 

4. All control surface hinge lines and inboard side edges are assumed to be 
represented by a sealed gap condition. 


5. -ATI control surface side edges are assumed to be parallel to the remote 

stream direction. 

6. Spanwise elastic twist distribution of the control surface hinge line is 
approximated by a cubic polynomial to allow a spanwise variation of the 
hinge line rotation © H (n). 

7. Control surface camber bending is treated within the regular lifting 
surface solution. 

8. Lifting surface solutions of planforms not having control surfaces may 
be obtained upon user option. 

9. Spanwise symmetrical or antisynmetrical analysis may be accomplished upon 
user option. 

10. Planform definition is assumed to consist of linear segments defined by 
input data. (See Appendix C on effect of locating downwash chords near 
junctions of linear segments of planform definition.) 

11. Downwash boundary conditions may be modified to include local streamwise 
velocity variations due to airfoil thickness effects. (See Appendix C.) 
(See reference 4 for limitations.) 

12. Input mode shapes may be defined at many points on the surface (such as 
obtained from a finite element structural analysis). 

13. Multiple k value and Mach number analyses may be made in a single 
machine run (see reference 4, section 3.7.2 LIMITATIONS). 

14. The printed output of the program, obtained on user option, consists of 
a definition of kinematic downwash matrix, mathematically discontinuous 
downwash matrix, residual downwash matrix, pressure series coefficients, 
chordwise array of pressures at prescribed spanwise stations, sectional 
generalized forces, and total generalized forces. 

15. Intermediate and/or final results may be optionally saved on tapes for 
economy of reuse. 
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VALIDATION STUDIES OF ASYMPTOTIC PRESSURE EXPRESSIONS 


Numerical investigations have been conducted to determine the validity 
of using any of the various asymptotic pressure expressions that have been 
proposed 3 ’ 7 * 8 for analysis of configurations having swept hingeline trailing 
edge controls. 

The criteria used to judge the correctness of the asymptotic pressure 
expressions is that the residual downwashes are required to be smooth and 
continuous over the lifting surface including the region of the control sur- 
face edges. Since the asymptotic pressure expressions are developed to satisfy 
the basic flow equation and the change in boundary conditions across the edges 
of the control surface, then the subtraction of the kinematic downwash dis- 
tribution from the distribution having identical discontinuity values must 
result in a residual distribution that is smooth and continuous. Otherwise, 
it must be presumed that the expressions are functionally insufficient to 
satisfy the original boundary value problem. 

Initial studies to determine the characteristics of the residual down- 
washes were conducted on a zero sweep planform having a 30% chord control sur- 
face. The analysis was conducted for several spanwise length control sur- 
faces. The pressure expression used in the analysis is the expression de- 
veloped by Landahl 3 for the zero sweep rectangular control surface and for the 
steady flow case of k=0, the expression is given as 

c p = - ^ to 9 [/ U-5 C ) 2 + e 2 (n-y s ) 2 - e(n-y s )] < 78 ) 


Analysis results obtained for all of the control surface configurations 
indicate that the residual downwashes are smooth and continuous over the en- 
tire planform and the final solutions provide generalized force values that 
are almost invariant with the number of collocation stations used in the 
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analyses. Chordwise pressure distributions obtained for spanwise stations 
located at large distances from the control surface side edges decrease in 
value montonically from the leading edge to the trailing edge without display- 
ing any chordwise waviness in any of the study configurations including solu- 
tions for very small span control surface configurations. 

It appears that solutions obtained for the zero sweep control surface 
case using the pressure expression of equation (78) will be almost variant 
with the number of collocation stations used in the solution process. 

Solution sensitivities were then obtained for the swept hingeline case 
using another loading function as suggested by Ashley 11 and as applied in 
computer programs developed by Darras and Dat 12 and by Zwaan 13 . 

The form of the k=0 asymptotic pressure expression is given as 


C p * - tog[/ (£-5 c ) 2 + e 2 (n-y s ) 2 - e(n-y s )] (79) 


where 


p * / ( l-M 2 cos 2 A )/ cos 2 A * / B 2 + tan 2 A 

C C v 


It should be noted that the only difference between this expression for swept 
hingelines and Landahl's expression for non-swept hingelines is the factor 3 
as contained in the multiplying coefficient. Unsteady lifting surface pressure 
distributions were obtained using the constant chord 25° sweep planform of 
reference 17 shown in figure 20 for which experimental lifting surface data 
are available. Motion of the configuration is composed of an oscillating out- 
board flap with all other components maintained in a stationary position. 
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Figure 20 - Experimental Planformof Reference 17 Showing the Slde-by-Side 
Control Surface Arrangement, (dimensions in meters ) 


Figure 21 represents the chordwise variation of the in-phase pressure 
component (ACpreal) for a spanwise station that is located at a large dis- 
tance away from the control surface side edge. Experimental values are 
represented by the open symbols and two theoretical curves are given that 
represent the solution results using two distinctly different distributions 
of collocation stations. The theoretical results for the (4,6) solution repre- 
sent the results obtained using 4 collocation stations on 6 downwash chords. 

The downwash chords are located at relatively large distances away from 
the control surface side edge. Theoretical prediction using the (4,6) dis- 
tribution follows the trends of the experimental values, except in the region 
of the hingeline gap where the gap is affecting the experimental results. 
Whereas, the solution results of (7,9) distribution contain a chordwise 
waviness near the trailing edge that is thought to be due to having 
the downwash chords located very near the control surface side edge. The 
presence of small pressure wave lengths of the (7,9) distribution results, in 
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Figure 21 - Comparison of Theoretical and Experimental In-phase Pressures. 


comparison with the (4,6) distribution results, indicates that the residual 
downwash distribution may not be smooth and continuous across the control sur- 
face side edge with the result that some of the low wave length assumed pres- 
sure modes are required to satisfy irregular boundary conditions being 
generated by the asymptotic pressure function of equation (79). 

Numerical investigation has been conducted to determine the character- 
istics of the residual downwashes being generated using the pressure expres- 
sion of equation (79) applied to the analysis planform of figure 22. Effect 
of the large hingeline sweep angle on residual downwashes should be well 
identified in using this configuration since tanA becomes equal to the same 
order of magnitude as e 2 . 

The spanwise variations of the residual downwashes obtained from the 
investigation for k=0 and M=0 are shown in figure 23. The spanwise 
variation of residual downwash for stations ahead of the hingeline is smooth 




( 100 , 100 ) 


( 150 , 100 ) 


Figure 22 . Planform Used to Evaluate Residual Downwash Characteristics. 

and continuous. However, for stations aft of the hingeline, there is a finite 
discontinuity as well as a slope discontinuity at the side edges of the con- 
trol surface. This side edge discontinuity may be analytically predicted, 
since it can be shown that the pressure expression will satisfy the boundary 
condition discontinuity across the hingeline, but will not satisfy the side 
edge discontinuity condition. Consequently, use of the pressure expression 
of equation (79) should be limited to configurations having zero sweep angles 
and large span length control surfaces with the restriction that the colloca- 
tion stations are spaced well away from the control surface side edges. 
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Figure 23 - Spanwise Residual Downwash Distributions Along Spars Forward 
and Aft of Hingeline Using Pressure Expression of Equation 79. 

Residual downwashes have also been obtained using the configuration of 
figure 22 and applying another pressure expression that has been used in 
computer programs developed by Rowe 1 *} Medan 15 , and Hewitt 1 ! The pressure 
expression contained in these references for k=0 is 

c p - - -j* &> g [ v < 5-? c ) 2 + BMn-yj) - 6(n-y s )] ( 8 °> 

This expression differs from equation ( 79 ) by having F included within 
the logarithmic function as well as appearing in the multiplying coefficient. 

The residual downwashes shown in figure 24 indicate that for chordwise 
stations ahead of the hingeline the residual distribution are smooth and 
continuous across the side edge. However, the distributions have a slope 
discontinuity across the side edges for stations aft of the hingeline. 
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It was also noted that for small control surface span lengths the magnitude 
of the distributions remained the same and the discontinuity peaks were spaced 
closer together which produces a highly localized oscillatory distribution in 
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Figure 24 - Spanwise Residual Downwash Distributions Along Spars Forward 
and Aft of Hingeline Using Pressure Expression of Equation 80. 


the residual downwashes. This localized distribution in residual downwashes 
is identified as being the major contributing factor that may produce 
questionable results for small span length control surface analyses. 

In order to alleviate the solution sensitivities caused by irregular 
residual downwashes new pressure expressions were obtained by reformulating 
the asymptotic expansion process following the procedures suggested by 

3 

Landahl. The step-by-step procedure is presented in appendix 'A'. The k=0 
pressure expression obtained from the new formulation is given as 

C p = - — log^j x 2 + B 2 (n-y s ) 2 - [e 2 (n-y s ) + x tanA]/^J (81) 

where X = S" x cs 


and the resulting residual downwash distributions are shown in figure 25 for 
the analysis configuration of figure 22. It should be noted that the residu- 
al distributions are now smooth and continuous for all stations either ahead 
or behind the hingeline. The chordwise pressure distributions at large span- 
wise distances away from the control surface were found to be montonically 
decreasing in value from the leading edge to the trailing edge. Furthermore, 
the final lifting surface pressure distributions remain almost invariant with 
the distribution and number of collocation stations used in the analyses. 
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Figure 25. Spanwise Variation of Residual Downwash Distribution Using 
Pressure Expression of Equation 81. 

Numerical investigations have also been conducted to determine the valid- 
ity of the equations developed to satisfy the change in boundary conditions 
across the side edge of an oscillating full chord control surface. 

Loading functions developed for the full chord control surface given by 
equation (65) have been used to evaluate the residual downwash distribution. 
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The analysis configuration is a constant chord 45° swept planform having 
a full chord control surface with side edges located at 60% and 80% of the 
semi span length. The residual downwash distribution is shown in figure 26 
and indicates that the residual downwashes are smooth and continuous across 
the side edges. 




y/s - Spamrise Station 


Figure 26. Spanwise Residual Downwash Distributions for a Full 
Chord Control Surface Hinged at 30 Percent Chord. 


It appears that the asymptotic pressure expressions developed for the 
swept hinge line trailing edge control surface and for the swept leading edge 
control surface does provide the proper downwash discontinuities and may be 
combined to provide analysis capability for the partial chord leading edge 
control surface. 



RESULTS, COMPARISONS AND DISCUSSION 


This section contains an illustration of the downwash subtraction pro- 
cedure and also provides pressure distributions and generalized forces result 
ing from analysis of five wing-control surface configurations for which ex- 
perimental data are available. The experimental configurations consist of 

(1) the swept-wing with a non-oscillating full span flap of reference 18, 

(2) a swept-wing having an inboard partial span control surface of reference 
19 for which non-pscillatory pressures are available, (3) a rectangular wing 
having an oscillating full span flap (reference 20), (4) a swept wing having 
oscillating side-by-side control surfaces (reference 17), and (5), a highly 

’ swept delta wing having leading and trailing edge controls (reference 23). 

Description of Downwash Subtraction Process 

The configuration of reference 17, shown in figure 27, is used to 
illustrate the downwash subtraction process. 



Figure 27.- Experimental Planform of Reference 17 Showing the 

Side-by-Side Control Surface Arrangement (dimensions 
in meters) 









The Inboard control surface is assumed to be oscillating about its hingeline 
at a reduced frequency of k = 0.372 and at a Mach number of H * 0. The rota- 
tion angle is assumed to be constant across the span of the control surface 
having a value of one radian measued in a plane perpendicular to the y axis. 
The remaining portion of the lifting surface is assumed to be at rest. 

The in-phase part of the kinematic downwash distribution obtained from 
the motion definition is shown in figure 28 for nine chords equally spaced 
across the semispan. The downwash has a uniform value of ^ = 1.0 over the 
control surface and is zero everywhere else. 



Figure 28. Kinematic Downwash Distribution Derived from 
Motions of Inboard Flap 

Downwash distribution resulting from the integral evaluation using the 
loading functions defined by equation (56) is shown in figure 29. The distri- 
butions are smooth in regions outside of the control surface and a unit change 
in downwash across the hingeline is obtained for the chords that lie on the 
control surface. — -- 
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i Since the modified distribution is smooth and does not contain any discontinu- 
ities, it appears that the pressure loading functions (equation (56)) that 
were developed from the asymptotic expansion process of reference 3 provide 
the proper description of surface loadings to be used in obtaining solutions 
for configurations having discontinuous downwash distributions. 

Steady-State Results for Full -Span Flap Configuration 

The full-span flap configuration of reference 18 is shown In figure 31 
for which experimental pressures were obtained for various combinations of 
flap deflections and wing angles of attack. The flap deflection and wing 
angle of attack were maintained at constant values for each experimental run. 



Experimental pressures were obtained along a streamwlse section located 
at the 50% semispan station. The longitudinal junction between the wing and 
flap was sealed to prevent flow leakage between the lower and upper surfaces 
at the hingeline. 
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The theoretical pressure distributions are obtained using a modi- 
fication on the boundary conditions (described in Appendix C) that in- 
clude local streamwise velocity variations due to airfoil thickness 
effects. A comparison of the experimental and theoretical results are 
shown in figure 32. 



TIP 

Figure 32.- Theoretical and Experimental Chordwise Pressure Distributions for a 
Full-Span Flap Deflected by 6 = 10°, a = 0.0° at M = 0.21 and k = 0 

The experimental sealed gap condition at the hingeline satisfies 
the theoretical assumptions (reference 3) and provides a one-to-one 
basis for evaluating the accuracy of the theoretical prediction method. 

The comparison indicates that the experimental values are theo- 
retically predicted within very close tolerances over the entire length 
of the chord even in the vicinity of the hingeline. 
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Steady-State Results for a Partial -Span Flap Configuration 


The partial-span control surface configuration shown in figure 33 
represents the experimental planform of reference 19 to obtain chord- 
wise pressure distributions due to a steady flap deflection. Pressures 
were obtained along a streamwise chord located at the 46% semi span 
station. The hingeline gap was sealed providing a one-to-one basis for 
comparing theoretical and experimental results. 



Figure 33.- Experimental Planform of a Partial -Span Flap Configuration of 
NACA RM L53C23 (dimensions in cm) 


Theoretical pressures were obtained along various streamwise chords 
spaced over the semi span and are shown in figure 34. The comparison 
indicates that the experimental pressures are accurately predicted by 
the theoretical technique over the length of the chordwise strip for- 
ward of the hingeline. The theoretical distribution over the control 
surface are only slightly larger than the experimental values. Con- 
sequently, it appears that the lifts and hinge moments may be predicted 
within reasonable accuracy limits for configurations having a sealed 
gap at the wing-control surface junction. 




ROOT 

Figure 34. Theoretical and Experimental Pressure Distributions for Partial- 
Span Flap Deflected by 6 = 10°, a = 0°, at M = 0.60 and k * 0 

Rectangular Planforms Having Full -Span Control Surfaces 

Two rectangular planforms having aspect ratios of 1.0 and 2.0 were tested 
in combination with a 40% chord flap for various values of reduced frequencies 
to obtain hinge moments and lifts due to the control surface motions. The 
experimental data are reported in reference 20. Each planform had a small gap 
along its hingeline that may have a small influence on the pressure distri- 
butions in regions near the hingeline. 

Experimental and theoretical values of hinge moments and phase angles 
as a function of reduced frequency are presented in figure 35 for the aspect 
ratio 2.0 planform. The theoretical hinge moments appear to agree with the 
experimental values; however, there is a discrepancy in the predicted phase 
angles that may be due to the open gap at the hingeline. 
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Figure 35.- Theoretical and Experimental Hinge Moments and Phase Angles for an 
Oscillating Full -Span Flap AR = 2.0 Configuration of WADC TR 53-64 


Comparison of wing lifts and phase angles are shown in figure 36. 
There is good agreement between the theoretical and experimental results 
over the entire range of reduced frequencies used in the tests. Conse- 
quently, it appears that wing lift is not greatly affected by open gaps 
at the hingeline. 




Figure 36.- Theoretical and Experimental Lifts and Phase Angles for an 

Oscillating Full -Span Flap AR * 2.0 Configuration of WADC TR 53-64 
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Experimental and theoretical results for the full span control 
surface configuration of aspect ratio 1.0 are shown in figure 37 and 
figure 38. Again, the theoretical method appears to predict reasonable 
values of hinge moments and total lifts; however, the predicted hinge 
moment phase angles deviate from the experimental values and the dis- 
crepancy may be due to the open gap at the hingeline. 




Figure 37.- Theoretical and Experimental Lifts and Phase Angles for an 

Oscillating Full -Span Flap AR = 1.0 Configuration of WADC TR 53-64 




Figure 38.- Theoretical and Experimental Hinge Moments and Phase Angles for an 
Oscillating Full -Span Flap AR = 1.0 Configuration of WADC TR 53-64 



The difference may also be caused by trailing edge flow separation 
that could have been experienced during the test since there is a very 
large adverse pressure gradient at the trailing edge that might not 
have been maintained during the experiment, (see page 315 of reference 
24 for the pressure gradient plot of the NACA 0010 airfoil section used 
in the experiment.) 

Effect of Hingeline Gaps on Chordwise Loadings 

A theoretical study has been made (reference 21) to investigate the 
effect of hingeline gaps on chordwise pressure distribution, and the 
results are shown in figure 39. 




Figure 39.- Non-Dimensional Lift Distribution (reproduced from reference 21) 



The study configuration is composed of a pair of two-dimensional equal 
chord lifting surfaces that are placed in a coplanar arrangement. The 
leading surface is at a zero angle of attack and the trailing surface is 
given an angle of attack by rotating the surface about its leading edge. 
Chordwise pressures are obtained for various separation spacings in steady 
flow. 


For a zero gap spacing, the loadings produced at the junction of the two 
surfaces exhibit the typical symmetrical logarithmic singularity identified 
for a discontinuous downwash distributions of a two-dimensional wing-control 
surface combination. As separation begins, the loadings on the wing near the 
wing trailing edge change from a singularity characteristic to a loading that 
becomes equal to zero. However, the loadings near the leading edge of the 
control surface change from a logarithmic characteristic to that of having an 
inverse square root singularity that causes a higher loading to exist on the 
control surface for the open gap condition than for the zero gap condition. 

Large separation distances between the surfaces result in loadings that 
are vastly different from the zero separation case. However, small separa- 
tions affect the pressure distributions only in localized regions near the 
hingeline. It appears that the zero gap theoretical method will provide 
reasonable predictions of total lifts and moments for configurations having 
small hingeline gaps. 

Side-by-Side Control Surface Configuration 

The side-by-side control surface configuration is shown in figure 27 for 
which unsteady pressures are obtained for various combinations of flap de- 
flections. The experimental model had small open gaps at the hingelines and 
at the side edges. Reference 17 provides no information to define the exact 
dimensions between the control surface side edges and the pressure measuring 
stations located in the vicinity of the side edges. The spanwise location of 
the experimental pressure chords were established by measuring the pressure 
chord locations off the planform drawing. 
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Figures 40 and 41 present comparisons of theoretical and experimental 
results for a mode shape in which both flaps are oscillating with the same 
phase and amplitude and the wing is maintained in a stationary position. 
There is good agreement between the theoretical and experimental results in 
all areas except in localized regions near the hingelines. The experimental 
values are less than the theoretical values forward of the hingeline and 
slightly greater than the predicted values aft of the hingeline. The varia- 
tions near the hingeline are attributed to the open gaps at the hingelines. 
Even though the comparisons deteriorate in regions near the hingeline, it 
appears that the total lifts and moments may be estimated within reasonable 
bounds for configurations that have small open gaps at the hingelines. 




Figure 40. In-Phase Pressure Distributions for Both Flaps Oscillating with 
the Same Phase and Amplitude. A. f = 0.82°, A - = 0.82°, 

A^ = 0°, k = 0.372, and M « 0 1,T * °* T - 
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Figure 41. Out-of-Phase Pressure Distributions for Both Flaps Oscillating 


Figures 42 and 43 present comparisons of theoretical and experimental 
pressure distributions for a mode shape in which the outer flap is oscillating 
and the inner flap and wing are maintained in a stationary position. 



HOOT 


Figure 42. In-Phase Pressure Distributions Resulting from Motions of Outer 
Flap. A. . = 0°, A„ - = 0.66°, A , = 0°, k = 0.372, and M • 0 
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Figure 43.- Out-of-Phase Chordwise Pressure Distributions Resulting from Motions 
of Outer Flap. A, - - 0°, A - = 0.66°, A = 0°, k = 0.372, and HaO 

i • T • 0 . T • W 


Again, the only area where the comparisons deteriorate is within 
the hingeline region and is probably due to open gap effects at the 
hingeline. 


Swept Delta Wing with Leading and Trailing Edge Controls 


The planform configuration shown in figure 44 represents the experimental 
planform of reference 22 where experimental studies were accomplished to 
investigate use of active controls to suppress flutter. Steady state hinge 
moments obtained within this investigation were reported within reference 23 
and are reproduced in figure 45. 
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Figure 44- Experimental Delta-wing Model of NASA TM X-2909 
(all linear dimensions in meters ) 
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Figure 45. Comparison of Measured and Calculated Static Hinge 
Moment Coefficients 


Two sets of analysis results obtained from the present procedure are 
superimposed on the test results of figure 45. The set of results denoted by 
the symbol ♦ represent the solutions obtained using standard flat plate 
boundary conditions. The set of results denoted by the symbol ^ represent 
the solutions obtained using local linearization of boundary conditions to 
include effects of incremental velocities due to thickness distribution as 
described in Appendix C. 

It appears that the application of local 1 inearization of boundary con- 
ditions does affect the hinge-moments of the trailing edge control surface. 
However, there is only a slight change in hinge-moments of the leading edge 
control surface due to local linearization effects. 

The analytical results appear to correlate well with the experimental 
data for the Mach range of interest. The comparison tends to deteriorate as 
the Mach number approaches 1.0 which may be due to shock effects that may be 
caused by the control surface deflection. 
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CONCLUSION 


’ A theoretical analysis and computer program has been developed for the . 

r prediction of unsteady lifting surface loadings caused by motions of leading 
edge and trailing edge control surfaces having sealed gaps. The program has 
been developed around a- systematic solution process where the discontinuities 
in the downwash distributions are separated (and handled separately) prior to 
applying standard lifting-surface solution techniques. 

Theoretical results are presented to demonstrate the versatility of the 
program capabilities. Comparisons of theoretical and experimental data are 
presented for five wing-control surface configurations. Two of the experi- 
mental configurations had sealed gaps at the hingelines, and two configura- 
tions had small open gaps. 

The comparisons indicate that reasonable pressure distributions are ob- 
tained for the sealed gap cases, but there are moderate differences near the 
hingelines for the open gap cases. 

Pressure distributions and hinge moment coefficients obtained for all 
trailing edge control surface configurations appear to approximate the experi- 
mental values within close tolerances. 

Hinge-moment coefficients obtained for the one leading edge control sur- 
face configuration does not appear to be sensitive to the use of local 
linearization of boundary conditions. Analysis results correlate well with 
experimental data for most of the Mach number range except for Mach numbers 
in the close vicinity of 1.0. 




APPENDIX A 


DEVELOPMENT OF PRESSURE EXPRESSIONS THAT SATISFY THE BOUNDARY CONDITIONS ON A 
TRAILING EDGE CONTROL SURFACE HAVING A SWEPT HINGELINE 


Pressure expressions are formulated such that the boundary conditions are 
exactly satisfied and the residual downwashes that remain after subtracting the 
downwash discontinuities are smooth and continuous across the side edge of the 
control surface. The analytical procedure used to obtain the expressions is 
the asymptotic expansion process suggested by Landahl In reference 3 . 

Solutions of the boundary value problem obtained by application of Green's 
theorem provides the required change in boundary conditions across the hinge- 
line and side edges such that subtraction of the discontinuous kinematic down- 
wash distribution result in residual distributions that are smooth and 
continuous for which standard lifting surface solutions may be readily applied. 

The analysis coordinate system is shown in figure 46 where the x,y 
coordinates are referenced to the inboard side edge of the hingeline. 


WING 



Figure 46 : Analysis Coordinate System 
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The Linearized boundary value problem is developed in terms of the 
pressure perturbation Cp = F e where the pressure magnitude satisfies 
the equation of flow 

e 2p_ + F~~ + P~ - 2ikM 2 P~ + k 2 M 2 F =0 (Al) 

xx yy zz x 


and boundary conditions on the surface of 

- 20 h [Mx-x c ) + 2ik - k 2 (x-x c ) U(y)] for x > x £ 

V , - 

0 for x < x c 

Scaling the x coordinate to remove the B 2 factor we let 



9 1 9 - 9 2 = l_iL_ 

3x 3 3x Q * 9x2 3 2 9 x q 2 


(A2) 


(A3) 


and the equation of flow becomes 


Vo + Vo + Pz o z o ' 2 0 ' x 


ikM 2 y + k 2 M 2p = o 


(A4) 


The terms of the boundary condition are modified by letting 


x-x c = x-y tan A 


= 3(x 0 -y 0 tan A/ 3) 
= e ( x o‘ y o tan V 


and noting that 


S(ax) = i <5 ( x) 
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Thus the boundary conditions are given as 

- 20 h [1 6(x 0 -y Q tan A-j ) + 2ik - k 2 e(x Q -y 0 tan A-| ) J 
0 


V 


x o t *c 


for x„ < x. 
o c. 


(A5) 


The first derivative term of (A4) is now removed by defining a variable 
ip def i ned as 


/ ikM 2 x. \ 

7 - exp (-g-Zj t 

which will change the flow equation into the form 

» + * . o 


(A6) 


Vo Vo Vo 


3 5 


(A7) 


and provide a boundary condition definition of 

- 20 h [^( x 0 -y 0 tanA 1 ) + 2ik - k 2 6 2 (x 0 -y () tanA 1 )J e -lkM2x o /e y Q > 0 


_ 

32 o 


y 0 < o 


To study the lifting surface problem in the vicinity of the corners all 
coordinates are scaled by e 


(A8) 


- x o - *o - z o 

x = - ; y - - ; * = - 


that modifies the differential equation into the form of 


\b-~ + il;-- + ill-- + £ : 

r xx Y yy v zz 


2M2 


k z M 


ip = o 


(A9) 


and provides a boundary condition definition of 

2 0 H (x-y tanA-j ) + 2ik - k 2 eB(x-ytanA-| ) J e” lkM eX ^ y — 0 




0 


y < 0 (A10) 


We now seek solutions of the form 


= t|> 0 e° + e + ^ e 2 + • • • • (All) 

and expand the exponential in the boundary condition and set to zero the 
collected terms of equal powers in e. We finally arrive at the following 
boundary value problems whose sum of solutions will satisfy the original 
boundary value problem. 

|”v 2 <l > 0 - 0 



‘ 20 H 
0 



x - y tan 



/ 


y > 0 
y < 0 


(A12) 


V 2 ^ = 0 


3 ^ 

3¥ 


- 20 h [2 ik - - k 6 (x - y tan A 1 )] y > 0 
0 - y < 0 


(A13) 


V 2 h - - 


k 2 M 2 
B 2 


4> r 



- 20 h [ 2 — g - x - k 2 6(x-ytanA 1 ) - k 2 M gi 6(x-ytanA 1 )j y >0 
0 y < 0 


L (A14) 

Since the boundary value problems contain (or can be made to contain) a 
Laplace equation V 2 \p=0 then use of Green's theorem may be readily applied 
to obtain solutions of the form 
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*<x.y,z> * - | J / / 


9^(x,y,0) 

3z 


dxd* 

((x-x) 2 + (y-y) 2 + z 2 y 


'z=0 


(A15) 


where 

30(7 J,0) 

3x 7=0 

represents the surface boundary conditions on the z»0 plane and the Green's 
function (developed from the method of images) 

v = 1 + 1 

((x-x) 2 + (y-y) 2 + (z+T) 2 f ((x-x) 2 + (y-y) 2 ♦ (z-7) 2 )* 

v satisfies the conditions of 
^ = 0 on 7 = 0 

dZ 

V 2 v = - 4ir6(x-x,y-y ,z-7) 

Integration of (A15) is accomplished in rectangular coordinate system 
that have coordinates defined within figure , and the coordinate subscripts 
and bars are removed with the understanding that these are scaled coordinates 
and that the solution will be finalized in physical coordinates at the end of 
this section. 


Solution of the Zeroth Order Problem 

The zeroth order boundary value problem is given as 
P.D.E. 

B.C. 


v\ = 0 

9'l' o (x,y,0) 20m 

= " 3 tan V 
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Substituting the boundary conditions into equation (A15) gives 


®H T T SU-y'tanA^Ufy^dx'dy' 


n f f 1 

*° (X ' J ' ,2) J J (^7)* +( x-r)H 2 ^ 


-OO -00 


Changing the lower integration limit of y" (due to U(y")) 
the x' integration produces the expression 


^ 0 (x,y,z) 


^ 00 
£h f 
Tre J 


_dy: 


((x-y^tanA^ ) 2 + ( y-y ") 2 + z*f 2 


substituting x=x-x c +ytanA^ in (A17) yields 


0 H r 

= is J 


<&1 


(((x-x c ) + (y-y") tanA^ 2 + (y-y") 2 + z 2 )^ 


Expanding the terms in the radical provides 


^ 0 (x,y,z) 


= ^ f 

ire J 




((l+tan^ ) (y"-y) 2 -2(x-x c )tanA 1 (y"-y)+(x-x c ) 2 


that is integrated to yield the function 


<P 0 (x,y,z) 


Q: 


H 1 


= + -s-r: 

irB ( a )' 5 


£(a)^ (an 2 + 2bn + c)^ + an + b j 


I -y 


(Ai6) 

and 


(A17) 


+z 2 ) % 


^(Regular) 

(A18) 
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where 


a * (1 + tan 2 A^) 
b = - (x-x c ) tan A-j 
c = (x-x c ) 2 + z 2 

The divergent upper limit is handled in the manner of Landhal and the 
expression is reduced to the form of 

<l> 0 (x,y,z) = Zn [(a) % (x 2 +y 2 +z 2 )^ - (y+xtanA,)] + ^(Regular) (A19) 

(a ) 4 


where x,y,z are scaled coordinates. A check may now be made on whether the 
boundary conditions are satisfied. The boundary conditions on the z=0 plane 
are evaluated as follows 


3 *0 _ ~ S H 1 

3z - -IS (J) % 


(a)'* (x 2 +y 2 +z 2 )"* + (y+xtanA^) 
a(x 2 +y 2 +z 2 ) - (y 2 +2xytanA 1 +x 2 tan 2 A-j ) 



(x 2 +y 2 +z 2 )^ 


reorganizing this expression by inserting a=(l+tan 2 A^) into the denominator 
of the bracketed term we obtain 


»♦. _ ~ e H 1 

U 716 (a) % 

and noting that 


j. (y+xtanA. ) 

(a) 2 + 

(x 2 +y 2 +z 2 ) 2 


lim 3^+TTT " *«(x) 
x+o 


(a) % z 

(x-ytanA-j ) 2 +az 2 
then 



Z-+-0 


" 0 H 1 T, vJ' y(l+tan 2 A,) 

778 (a) % _ 3 + |y|((l+tan 2 A 1 ))^ 


ir6(x-ytanA 1 ) 


3*0 

W 


(A20) 


= - 2©j_| 6(x-ytanA-j) U(y) 

z->o — 

which is identical to the scaled boundary conditions of equation (12-A) and 
thus the function of equation (A19) is a valid solution of the zeroth order 
boundary value problem. 

Scaled coordinates are used in equation (A19) and the expression may be 
returned to physical coordinates by replacing 

x by y by and z by | 

tanA-j = ^ -| ; B 2 = B 2 + tan 2 A 

to provide the zeroth order pressure expression on z=o of 



(B 2 y+xtanA) 

6 


+ * 0 (Regular) 


(A21) 


Solution of First Order Problem 

The first order boundary value problem in scaled coordinates is defined 
by 

[ V 2 ^ = 0 

-20 h ^2i k - i kM 2 x 6(x-ytanA^)j y > 0 

0 y < 0 (A22) 


3 *! 

W 
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and is divided into two separate problems corresponding to the two terms within 
the boundary condition statement. The first part of the first order problem is 
defined as 

= 0 

= - 20 h [21k] U (x-ytan A ] ) U(y) 
z=0 


v \ 

8 \ 


Inserting the boundary conditions into equation (A15) and changing limits 
due to the unit functions provides the expression 


(x,y,z) 

'a * 


_ , , oo oo 

°H’ k // 


dx^dy ' 


o y'tanA-j ((x-x') 2 +(y-y') 2 +z 2 y 


(A23) 


A straight forward integration of this equation followed by a transformation 
of the coordinates into physical coordinates z=o yields the solution 


eih (x,y,z) 3 

'a 


0 H 2ik 


~ r,~> v /3 2 y+xtanA,V 

y £n[(x 2 +3 2 y 2 ) a - xj + —■ c i n l_j 


(A24) 


The second part of the first order problem is given as 


r V 2 ^ = 0 


3 * 1 , 

3z~~ 


20 u i kM 2 


H 


z->0 


p — 6(x-ytanA 1 ) U(y) 


Again inserting the boundary condition into equation (A15) and performing 
the indicated integration yields the expression for the pressure term in 
physical coordinates on z=o plane given as 




B 


0^1 kM 2 tanA 
ttB 2 B 3 


(xtanA+yB 2 ) In (x 2 +B 2 y 2 


_ / B 2 y+xtanA 

V 3 



(A25) 


The above solutions of the first order boundary value problem have been 
obtained using the same techniques of Landahl in the evaluation of the integral 
at the divergent upper limit. 

Solution of the Second Order Problem 


The definition of the second order boundary value problem given in scaled 
coordinates is 


= - 


2m2 


3 2 % 


3*2 

3z 


= 20 . 


Z=0 


KW^^-tanA^-Z^Utx-ytanA,) 


+ k z 3(x-ytanA^ ) U(x-ytanA^) 


u(y) 


(A26) 


The differential equation is a Poisson equation whose solution requires 
some added manipulations over the procedures used in obtaining the preceding 
boundary value problem solutions. The technique applied here to obtain solu- 
tions of the second order problem is one of redefining the boundary value 
problem into separate parts and obtaining individual solutions whose sum will 
satisfy the original boundary value problem. 


We first seek solutions that satisfy the following boundary value 
problems 


V 2 i|> 

3^2 


PL . 

2 

PL I 


9z 


z =0 


k 2 M 2 

3 2 


= 0 


*o 


(A27) 
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V 2 <|/ 2 A = 0 


W 


z=o 


fi k 2 M 4 x 2 

6(x-ytanAj ) U(y) 


(A28) 


V 2 ^ 6 = 0 


3^2 

W 


= - 40 h k -- ^ -- X - U(x-ytanA^ ) U(y) 


(A29) 


v 2 ^ 2 c = 0 


3^ 

W 


= 20m k 2 $(£-ytanA, ) U(x-ytanA-,) U(y) 
z=o 


(A30) 


Solutions of the last three boundary value problems may be obtained in the 
same manner that was used to obtain the solution of the zeroth order boundary 
value problem. Insertion of the boundary conditions into equation (A15) and 
performing the straight forward integration will result in the desired solutions 
which will be summarized at the end of this development. 


Solution of the Poisson equation and its associated boundary condition is 
obtained by seeking a solution of the form 


4>‘ 


, Pq * *, p ♦ *, q 


such that ij >2 is a solution of 


vy * - 
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p p 

Once Is obtained then the boundary conditions developed by i|* 2 are 

determined by evaluating the function lim ( am p / 3? ) 

z -*0 ' v 2 / / % 


Then a new boundary value problem is formulated such that ij> 2 q will satisfy 
V 2 ^ 2 q = 0 


3*0 


32 


= - l^im 
z-*-o z-*o 


3*2 

dz~ 


Thus the sum of ^ 2 and \p^ will then satisfy the differential equation 

V ■ - ^ *0 


having boundary conditions 


dip. 


Pq 


32 


= 0 


z-*-o 


In order to simplify the expressions a coordinate transformation is made 
such that new coordinates are defined along and perpendicular to the hinge 
line as shown in figure 47 . 



Figure 47 : Rigid Rotation Coordinate Transformation 



(A31) 


x = x cos A-| - y sin A-j 


x = x cos A^ + y sin A-j 

y = x sin A-j + y cos A^ 

or 

y = - x sin A-j + y cos A^ 

1 

N| 

11 

N? 

1 


z = z 


Using the above relationships it can be shown that 
V 2 (x,y,z) = V 2 (x,y,z) 


(x 2 +y 2 +z 2 ) % = (x 2 +y 2 +x*) % 


= R 


and that the solution of the zeroth order problem may be expressed as 


QmCOsA, 

^ 0 (x,y,z) C R ~£l 


(A32) 


We now seek a solution of Poisson's equation by developing an expression having 
all the required characteristics that remain after performing the three co- 
ordinate second derivative operation. One such expression that may fulfill 
these requirements is 


i|> 2 P = AR 2 In (R - y) + BR 2 (A33) 

where A and B are to be determined such that 
V 2 i|> 2 P = C Jbi (R - y) 

where c = 

now 

V 2 i|> 2 P = A V 2 (r 2 -6i(R - y)) + BV 2 R 2 (A34) 
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and noting that for a general function (uv) 


V 2 uv = uV 2 v + 2 (vu) • (Vu) + v • V 2 u 


then 

V 2 ^ 2 P = A j^R 2 V 2 £w(R-y) + 2(VR 2 ) • (v^w(R-y)) + fn(R-y)V 2 R 2 ] + BV 2 R 2 

now V z Jbi{ R-y)=0 since this is a solution to the zeroth order problem, also 
V 2 R 2 =6 and 2 (VR 2 ) • (v£n(R-y )) = 4. Therefore, the expression becomes 

V 2 ^ P = A [4 + 6 U (R-y) ] + 6B 
= 6A Zn (R-y) + (4A+6B) 

and constants must take on value of 



to satisfy the differential equation. Finally, the pressure expression that 
satisfies the Poisson equation is given as 



k 2 M 2 \ e H cosA l 
. 3 2 / ttB 


§1 b, (R-y) - 


(A35) 


or 


^ - 2 ] 


C 


k 2 M 2 \ Q H COsA 1 
. e 2 / tt$ 


p 

The boundary conditions generated by ik, are now determined by performing 

P c 

the differentiation of ij> 2 with respect to z and evaluating the limiting 
functional value as z -*• 0. 


dip. 


3Z 


■ Tb{ 2z [ 3 *’ (R ‘ y) ' 2 ] * 3r! [trVt t]] 

= — f< 

18 r 


6z£w(R-y) - 4z + |^iyy | 


(A36) 


The limiting value of the first two terms as z 0 are identically equal 
to zero and we are left with the limiting value of the last term. 

«.£-»■ KBH”" f&T - S (£=,><**> 


= gir6(x)(x 2 +y 2 )^ [(x 2 +y 2 ) ?s + y] 

r° ■ 


= §w6(x) |y| [ |y| + y ] 58 j 2 

l 26(x)|y| y > o 


y < 0 


r 

lim = £ir<S(x)y 2 U(y) 

I-o 8? 


(A37) 


In summary, the function 


*2 P = T§" ' 2 ] 


is a solution of 


V 2 * 2 P = an(R-y) 


and produces the boundary condition given as 

lim — — = c56(x)y*U(y) where C =1 

1-0 3? 

fk 2 M 2 \ °h cosA 1 
k 8 2 / tt3 



The next step Is to find a function such that ft satisfies the 

boundary value problem 

V 2 ^ = 0 




(A38) 


Following the same procedure that was used to solve the zeroth order prob- 

p 

lem results in the expression for i|> 2 given as 




1 ' 

I yft - [f 1 - 1 y 1 ] tnlR-y) ♦ ? 


Finally, the solution of the Poisson equation 
^ *0 

having zero boundary conditions 

3^21 


3z 


= 0 


z=o 


is given by the sum of solutions 
4> 2 = ♦; 2 + 


(A39) 


*2 s ? (R 2 -y 2 )MR-y) - §6 (4R 2 + 9yR + 6y*) 


(A40) 


where 



I 

The expression for ^ 1S separated into singular and regular expressions 
that are functions of the scaled and transformed coordinates 


*2 55 f(R 2 -y 2 )^( R -y) + ^(Regular) (A41) 

The expression given in physical coordinates is obtained by first rotating 
the coordinate system back into the scaled coordinate system (where 'o' sub- 
scripts are used to denote scaled coordinates) then transform from scaled 
coordinates to physical coordinates. 

R = (x 2 +y 2 +z 2 )^ = (x 0 2 +y 0 2 +z 0 2 )^ 

y = x Q sinA 1 +y 0 cosA 1 


The singular part given in scaled coordinates is 

^2 (x o’ y o’ z o ) 4( Cx o 2 % 2+z o 2] - [x o sinA l% cosA i ;|2 ) £n [ i:x o 2+y o 2+z o 2]%_Cx o s1nA l +y o cosA l 

(A42) 

The expression in physical coordinates is finally developed by applying the 
scaling transformations 


'0 = fe ; *0 " e ; z o = t E 0 : tanA l = 


tan A 


Thus, the singular part of the pressure expression obtained from the solution 
of the Poisson boundary value problem is given as 

e 2 iL(x,y,o) = — ~ (x-ytanA) 2 £n f(x 2 +3 2 y 2 ) % - / l!M ta n A\] (A43) 

* 4tt$ 2 8 3 L \ 3 /J 

The rest of the solutions required to satisfy the remaining second order 
boundary value problems defined by equations (A28) , (A29) , and (A30) may be 
obtained by Inserting the boundary conditions into equation (A15) and perform- 
ing the indicated integration. The solutions obtained and expressed In 
physical coordinates are summarized as follows: 


I 


V 2 * 2 A 


= 0 


dlp 2 

JF 


0 H k2M V 


«(vy 0 t a "V u(y o) 


V° 


BOUNDARY VALUE PROBLEM 
IN SCALED COORDINATES. 


Solution given in physical coordinates is 


A © H k 2 M 2 tan 2 A 

e 2 ty 9 = — 

4tt 3 2 B 5 


(x-ytanA) 2 - 2(^f^ + ye) 2 £n[R 2 -y 2 ] 


where 


g 2 +xtanA 


R ? = (x 2 +e 2 y 2 ) 2 ; y ? ■ 

L 3 


V 2 iJ; 9 B = 0 


9z" 


2 2 

• 40 hV*~ 


v° 


BOUNDARY VALUE PROBLEM 
IN SCALED COORDINATES. 


Solution in physical coordinates is 



20 H k 2 M 2 

Ttg 2 


(x-ytanA) 2 

26 3 


1 

3 


(x 2 -xytanA) 


MR 2 -y 2 ] 


20 u k 2 M 2 

ip — W -en[R 2 -x] 


(A44) 


(A45) 


(A46) 


(A47) 
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BOUNDARY VALUE PROBLEM 
IN SCALED COORDINATES. 


V 2 ^ 0 = 0 

= 20 H k 2 $(x o -y o tanA 1 )U(x o -y o tanA 1 )U(y o ) (A48) 

v° 


f- (£|^- xy)i„[R 2 -x] 

The pressure coefficient that satisfies the boundary value problem of 
equations (Al) , (A2) may now be expressed as a sum of the individual 

solutions 


Solution in physical coordinates 


0 

e 2 *? C = JL_ ( x-y tanA ) 2 £n[ R« -y 9 ] - 

* 2tt6 l l 



P(x,y,o) = (ip 0 + ei^i +e z ip2 ) e 


w 


(A49) 


F(x.y.o) = Cl-£n^(x 2 +B 2 y 2 )^ - ^^tanA^ j + C2 y £n j(x 2 +e 2 y 2 ) % - xj 


(A50) 


where Cl = 


- f l+2i 

irg l 


1 +2i k (x-y tanA) - — (xtanA+y 3 2 ) - — 1 — -( x-ytanA) 2 

B 2 B 2 4B 2 B 2 


k 2 M H tan 

4 B 2 6 

2k 2 M 2 
B 2 


|S^[(x-ytanA) ! -2P|^ + ey) a ] 

( x-ytanA) 2 -(x 2 -xy tan A)j - j- (x-ytanA) 2 j e 02 


C2 = 


. k ^ x + kH^yl e 




in 


The expression for F(x,y,o) may be simplified by expanding the 
exponential, carrying out the indicated multiplication and collecting terms 
that are linear in x and y. The resulting expression then becomes 

F(x,y,o) - - ^ Fl + 2fk(C-C c ) + : ^C-5 c )j in |(x 2 +e 2 y ! ) % - Bi ^_ ta - nA j 

- ^ [21k - k 2 (£-? c )] y 2»[(x 2 * e 2 y 2 ) % - x] (A51) 

The limiting value of this expression for A s 0 is equivalent to the 
expression developed by Landahl (reference 3 ) for the rectangular non-swept 
control surface. The above expression for A=0 is 

P(x,y,o) = - — 1 + 2ik(l--5 

t re L 

-Sl[21k-k 2 ({-c c ) ] y In ^((C-C c ) 2+ S 2 y 2 )^ - (?-5 c )j (*52) 

Thus, the expression for ^(x.y.o) appears to have the necessary 
characteristics to provide a proper solution of the swept hingeline boundary 
value problem and compares favorably with the widely accepted expression of 
Landahl for the zero sweep case. 


c ) ♦ in ((£-£ c ) 2 +S 2 y 2 y - Sy 
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APPENDIX B 


.DEVELOPMENT OF PRESSURE DISTRIBUTIONS TO SATISFY THE BOUNDARY 
CONDITIONS FOR A SWEPT LEADING EDGE CONTROL SURFACE 


Pressure distributions are formulated using the asymptotic expansion 
process suggested by Landahl (reference 3 ). The analysis coordinate system 
shown in figure 48 represents a segment of a swept wing leading edge having 
a control surface that oscillates about a hingeline aft of the leading edge. 



Figure 48: Analysis Coordinate System 

Solutions of the mixed boundary value problem are obtained using the 
linearized partial differential equation of flow given as 

3 2< |>-- + <P yy + <t> zz - 2ikM 2 <j>- + k 2 M 2 <j> = 0 (Bl) 

The motion of the control surface is assumed to be sinusoidal and is 
defined in non-dimensional coordinates as 

z c • 0 H (x - x c ) e 1kt 
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The boundary conditions resulting from the control surface motion are 

y<0 on z=0 

x>y tan A 

y<0 on z=0 (82) 

x>y tan A 

x<y tan A 

(ahead of leading edge) 

The first two expressions may be combined into a single expression for the 
change in boundary condition across the control surface side edge 

<t> z = (z + ikz ) U (y) 
x 

<J> z = © H C(l-ikx c ) + ikx] u (y) (B3) 

where U(y) = {J 


H*-^ +1kz c for 

H ■ 0 f ° r 

4. = 0 


A transformation or scaling of coordinates is made to simplify the differential 
equation by defining 


= x/B; 



z 


that results in the flow equation 


r x_x_. 
o o 


In the 


+ i y + 4> z z - + k 2 M 2 <J) = 0 

y o y o z o z o 13 *0 

new coordinate system the equation for the leading edge becomes 


(B4) 


or 


V = y 0 tan A 


y = w tan A 
o y o 3 
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Defining a new angle, A-j , such that tanA-j=^|^ then the leading edge 
equation becomes x 0 =y Q tanA 1 and the boundary conditions In the new coordinate 
system are 


<f» = 0 

for 

V y o tanA l on z ~ 0 

H ■ O+vx 0 )U(y 0 ) 

for 

x 0 >y Q tanA^ on z= o 

where X = 0^(1 -1kx c ) 



v = 0 H ike 




(B5) 


The first derivative term in equation (B4) is then removed by defining a new 
variable, 4>, defined as 


ikM 2 x /B 

= e 


(B6) 


Insertion of equation (B6) into (B4) results in the transformed differential 
equation 


’J'v 


+ ^ 


+ V 


z o z 


. * ("i 


ip = 0 


oo o' o 

The transformed boundary conditions are then defined as 


(B7) 


iji 5 0 


||- - (X+vx Q )e 


-ikM z x 0 /B 


U(y 0 ) 


x 0 <y 0 tanA 1 on z=d 
for * 0 >y o tanA l 


(B8) 


In order to study the problem in the vicinity of the leading edge corner all 
coordinates are stretched such that 


In terms of the stretched coordinates the boundary conditions become 

tj> = 0 x 0 <y Q tanA 1 on z Q =0 

„ , -ikM 2 ex /B , % 

||- * e(A+vex 0 )e U(y Q ) V y o tanA l on z Q =0 (B9) 

0 

Expanding the exponential and collecting terms of e n the second boundary 
condition is then defined as 

||- - eXU(y 0 ) + e’(v - x 0 U(y 0 ) + 0(e ! ) (BIO) 

The variable ifi is also expanded in a series given as 
ijj * + ^i e + ^2 e2 + ‘ ' (Bll) 

and inserted In the differential equation and boundary conditions terms of 
like powers in e are collected and define boundary value problems for which 
solutions are to be obtained whose sum will satisfy the flow conditions at 
the leading edge corner. 


V 2 ip 0 = o 

♦o = 0 



V 2 ^ = 0 


*1 


= 0 


3^ 

9z o 


xU(y 0 ) 


(zeroth order problem) 
on z Q =0 x Q <y 0 tanA 1 

on z Q =o x 0 >y Q tanA ] 

(first order problem) 
on z Q =o x Q <y 0 tanA 1 

on z Q =0 x 0 >y Q tanA 1 


(B12) 


(B13) 
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(second order problem) 



\p 2 = O on z=o; x Q <y 0 tanA 1 


= (V - x 0 U(y 0 ) on z=Oi x 0 >y„ tanA, 


( 814 ) 


Solution of the zeroth order boundary value problem results In an 
expression that is given by 


*o = C 


7——! A+ ^ (x-ytanA) 3 * 

\(x-ytanA)*/ 8 2 


(BIS) 


where the coefficient, C, can only be determined by a global integration 
over the planform. Since the function contains a square root and inverse 
square root term that are available in the main lifting surface solution then 
this zeroth order solution will be ignored. 


Solution of the First Order Boundary Value Problem 


A rigid rotation of coordinates is made such that the coordinates are 
parallel and perpendicular to the leading edge as shown in figure 49 . 



Figure 49: Rotated Coordinate System 


The new coordinates are defined by 


x = x Q cos A-j - y Q sin A 1 
y = x Q sin A ] + y Q cos A ] 

z = z Q (B16) 

and the old coordinates given in terms of the new coordinates are 
x„ = x cos A, + y sin A, 

0 1 ' (B17) 

y Q = - x sin A^ + y cos A^ 


The above transformation is a rigid rotation without stretching so that 
the orthogonality of the axes and scales are preserved. Thus the transformed 
first order boundary value problem is given as 


’h + ^1 + ^l =o 

xx yy zz 


ip-. = 0; x<o on z=o 


3^1 

gj- = Xllfy-xtanX^ ) ; x>o on z=0 


(B18) 


A Fourier transform in the y direction is applied to reduce the problem to a 
two dimensional problem with a parameter. The transformed differential 
equation becomes 

'f'-i + 4 ;, - s 2 iL = 0 (B19) 

xx zz 


118 


Transformation of the off wing boundary condition yields 
“ 0 x<o on z=o 



The second boundary condition is transformed by 


__ J e ->sy (J(y-xtanAj)dy 


- / e" 1sy dy - - 

J 4 


xtanA, 


is(2ir)' 


on z=o 


The convergence at the upper limit is circumvented by considering s to be 
complex having lmag(s)<0. Then the boundary condition is 


z i s (2tt)" 


e‘ isxtanA l ; x>0 on z=o 


(B20) 


Following Landahl in the construction of a solution of equation 
solution is assumed to be of the form 


a trial 


*1 ■ A / 


» - |s|r 1 -is(x-x')tanA 1 


1 1 


(B21) 


where © 1 = tan -1 and r-j 2 = x' + z 2 

The quantity A is determined by having ^ satisfy the boundary con- 
dition on the lifting surface (equation (B20) ). Thus, the above expression 
is differentiated with respect to z by using the expression for the derivative 
in polar form 


9 - C0SGl 9 + <inn 9 
iTj - 90 + sin0 l 
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x -|s|r.j-is(x-x')tanA.| 


3 *i a r "'''T'*' /w, “i f e, e, I 

3z — 7 J ziTt [ cos9 i s1n r + cos r s1ne i | i% ' 




- A l s l / 


x -|s|r-j-1s(x-x')tanA^ 


sin0-j cos^r- dx' 


(B22) 


-oo* 


The second integral vanishes for z-*o since 0^=0 or tt. The derivative 
becomes 


8 ^ 
3 T 


Z=0 



-isxtanA^ 



- |s |r^+isx"'tanA 


3/2 


30, 

sin-^- dx" 

4 


(B23) 


The apparent singularity is removed by adding and subtracting 
within the expression to obtain 


sin j 0 1 

7Jil 



Landahl has demonstrated that the second integral vanishes for z-»o. Since 
9-|-»o for x>o and 0 -|-*it for x<o, then the limit as z+o of the right 
hand side becomes 
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I 


3z 


z=o 


= + a 
+ 2 e 


-IsxtanA 


# r -|s| |x'|+1sx'tanA, t 

1 f Le 1 - lJ 

L ixi 3/2 

. -IsxtanA. « f -x'[|s|HstanA,] 1 

= + r e J vi ( 

a .. -» 


dx' 


Integrating this expression by parts yields 


dz 


z-o 




-isxtanA 


'{- 2 [e' xCi,MsUnA ' ] -,]. ( x'r 
00 -x c |s |+istanA 1 ] 


2 [ | s | +i stanA^ ] J 




dx" 


) 


Since 


w 

/ 


e" ax dx' 


0 X 




* 2 / .-»* dt - (0 


h 


then the expression for the derivative is 


3^ 


3z 


-isxtanA, 

= - A (tt[ | s | +i stanA, ] ) 2 e 
z=o 


This expression for the derivative must be equal to the boundary 
equation (B20) that results in evaluating A as 


A = - 


A isir(2[ |s 1+istanA-j])' 


(B25) 


(B26) 


(B27) 

condition 

(B28) 
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Finally i//, becomes 



-isxtanA, 

X e 1 

Itts (2C l s 1+istanA^ 3)^ 


/ 


-CO 


-Isjr^+isx'tanA^ 



cos 



(B29) 


Taking the inverse transform of equation (B29) and interchanging the order of 
integration yields 


*1 


2tt 


X 

m 


/ 


COS 


© 


°°-iy is[y+(x'-x)tanA,]-|s |r, 

f § , ds 

-»-i Y is[|s|+istanA 1 ]* 


dx' ( 630) 


Now, with p = y + (x ' - x) tan A-j , the inner integral becomes 


F(P>- / 


-iy isp- 1 s | r 1 


-oo-iy 


is[|s|+istanA^]" 


it; ds 


(B3I) 


The integrand contains a non-integrable singularity, however, this integral may 
be evaluated if the singularity is softened by first taking a derivative with 
respect to p 


<»-iy isp-|s|r-| 

F‘(P) = f £===^- ds (B32) 

_oii Y VM +istanA l 


Then F'(p) may be evaluated for the limiting case of lim F'(p) and 
obtain Y-K) 


oo 

F‘(P) = / 

— OO 


is -|s|r-j 
e 

^|s|+istanA^ 


ds 


(B33) 


where the path of path of integration is along the s(Real) axis. 





F' 


(p >■•/ 


isp-|s|r. 


o isp - 1 s r-, 


' -ys+istanA^ 


ds + 


f 


•^s+istanA^ 


ds 


or 


oo isp-sr, 0° -isp-sr-. 

F 1 (p) = / , 6 — ■■■ - ds + f == 

J Y s+1stanA i 0 V s_lstanA 1 


ds 


r e" Sri r e isp e " isp 

= J -=- y-f ■- + ? == ds 

0 V s |_vl + ' itanA i yl-itanA-j 


Observing that 1 +ltanA-j = \/ l+tan 2 A-j e^ = secA-je^ then upon 
expressions we obtain 


F * (p) = 2^cosA^j 


F' (p) = 2^cosA 1 
= 2VcosA 1 




(B34) 


combining the 


(B 35 ) 


(B36) 
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Now F(p) may be obtained by an Integration of F'(p) as follows: 



(B37) 



(B38) 


Since the lower limit, y, is an arbitrary limit then value of y is 
selected such that the second term vanishes; that is, we select y such 
that 



(B39) 


now 


r-j-iy = ^r 1 2 +y 2 




124 



Thus y is selected such that y = -r^ tan and consequently F(p) is 
given as 

( -iA-,/2 , ] 

e yj r^-ip J (B40) 

and the expression for ^ given in equation (B30) may be put in the form 

(B41) 


*1 ' - in 3 n 


f mirni F(p) dx' 

i Vn" 


or 


♦.-iv 


cos (0-j /2) r -iA,/2 
Im e 




dx'l 


and since Im(Z)= Im(x-iy) = - y = - Im(x+iy) then ifij becomes 


. _ 2 x r 

♦i - - — j 


2X ~ cos(0.|/2) ycosA^ 


f r T 


Im 


r iAn/2 ' 

e -yr^+ip 


dx' 
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now 


replacing p by p = y + (x'-x) tan A^ 


. _ 2x r 

“ ” u J 


cos(0-|/2) 


Im 


f I'V 2 , 

| e -yr^cosA-i 


+i (ycosA-|+(x'-x)sinA 


,)Jdx' 


Now let n = y cosA-j - x sinA^ and ^ becomes 


. _ 2A r 

*i - - r J 


2X * cos (0^/2) 


V 7 ? 


Im 


|.' v! V 


r^ cosA^ +1 x ' s 1 nAj +in 


in | dx' 
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The value of i|>, on z-*o may now be established by noting that lim r, = |x'| 

z->o 

and that the integral from to 0 is equal to zero since 0^ = ir for x'<o 
on z=o, thus the expression for on z=o becomes 




Yx'cosA-j+ix'sinA-j+iri jdx' 


(B44) 


Rather than attempting to evaluate this integral in its present form introduce 
a new function ^ such that 


+ 1 ^ ■ 


2 \ JV 2 

— e 


/ 


yjx 1 cosAj +1 x ‘ s i nA^ +i n 


dx' 


(B45) 


and note that ^ has the same definition as before in that it is equal to the 
imaginary part of the integral. 

iA,/2 i A-. -iA,/2 

Now with e - e . e 



(B46) 


The integral may be evaluated by making a substitution of x‘ = jj- and 
integrating the expression by parts to yield 


♦ 1*, - - a(e lAl V^i ♦ ) 

- iA l 

where a = ine and n = y cos A^ - x sin A^ 


(B47) 
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The function ^ is the Imaginary part of the equation and is given in 
stretched rotated coordinates and may be expressed in the original stretched 
coordinates by transforming back out of the rotated coordinate system by 
using 


x = x Q cos A-j - y Q sin A 1 
y Q * y cos A-j - x sin A 1 = n 


^ 



V ( x 0 cos A 1 -y-j sinA-j ) (x Q +iy 0 ) cosA 1 



Vv iy 0 + V x o~y 0 UnA i 

VV 1y o ‘ Vv y o tanA i 


(B48) 


Now making a new definition q = y Q /x 0 and replacing x Q by x/Be where x 
is physical coordinates, the above equation becomes 


+ 1^ = 


l 

+ ^ [ p pv 

d [yniq - 


(B49) 


where B-j = tan A-j 

Then may be obtained by evaluating the imaginary part of the above 

equation 


= * Im j- f f e 1 V( 1 -B 1 q)0+iq) 

[ V 1+i ‘q + V 1_3 i q l 1 

Vi+iq - V 1_3 i q Jj 


+ ^-lo 9 
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After some manlupulations the expression becomes 
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2A x 
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Secondary Order Boundary Value Problem 


The second order boundary value problem may be stated as 


V 2 *o = 


2 m 2 


k M 




* 2 = 0; z=o and x Q <y 0 tan 


3^2 

ST * c 2 x o u(y o ); z=0 > *0 y o ta " 1 
0 


where c 2 - (v - iJfx) 
v = i 0 ^kB 
A = © H (l - ikx Q ) 

This problem is separated into two parts by the following procedure 

~ z . _ k 2 M 2 . 
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(B51) 
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$ 2 i |> 2 = 0 

^2 = 0 on z=o; x Q <y 0 tan A-j 

3\pp 

= C 2 x Q U(y Q ) on z=0; x Q >y 0 tan A 1 (854) 

Since it is indicated in equation (B12) that \p Q contains functions that 
do not provide for any change in boundary conditions across the control surface 
side edge then it follows that is a regular function and thus the solution 
of the boundary value problem of equation (B53) must also be a regular function 
that will be obtained from the global solution portion of the solution process. 

A solution, Tiig » equation (B54) is taken to be given as 
C x o 

^2 = T~ J (n»y 0 »z 0 )dn (B55) 

6 l y o 

where B-j = tan 

The function ^ satisfies the differential equation since i|>j was obtained 
on that basis. Also, ^ 2 =o on z Q =o for x 0 <$-|y 0 since ^ = 0 for 

V 6 lV and 

3^2 

5^ 2 . 0 = C 2 <X o ' e l y o’ U(y o> 

due to the fact that 

3>^i 

* x U <*o> for x o > 6 1 » o ' 

0 Z-0 

The second order solution, may be obtained using a modified form of 

equation (B48) to define ^ within the integrand 
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To shorten the notation let 
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Replacing e 'cosA^ by and evaluating the two integrals results in 


the expression 
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Upon collecting like terms the expression becomes 


% = 


t 


'2 Tm I ? V y o (H36 l > 

IT Im 1 2(l-ie 1 ) v 


where 


v o a V x o + fy 0 
V 1 " V x o - Vo 


o v l + *o( <x o- (1 f’ 0 *o) l °s ( v^-)J 


(B59) 


J 



Thus, expressions for the solutions of the first and second order boundary 
value problems are now defined by equations (B48) and (B59) , and it only 
remains to obtain the pressure coefficient in terms of these solutions. By 
definition, the pressure coefficient C p is given in physical coordinates 
x, y by 

Cp(x,y,o) = - 2|<j>- (x,y,o) + ik<f» (x,y,o)j ( B 60) 

where <j>(x,y,o) obtained by the e expansion process is given as 
4 >(x,y,o) = e lkM2x/e2 (t|> 0 + zipf + + • • • ) 

= 2 |> MW (B61) 
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Now \J>i may be restated in stretched and scaled coordinates as being 
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and the expression given in physical coordinates (x,y) is obtained by 
replacing x Q by y Q by ^ and B-j by to yield 


= ‘ r- Ira |i 6+tanA v o 
where v = V x+i Sy 


,v?M^ 

V ~V, 
0 1 - 


(B65) 


v-j = yj x-ytanA 


Equation (B65) along with its derivative with respect to x may be inserted 
into equation (B63) to give the pressure coefficient due to the first order 
solution as 
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J m ([B+1tanA]v Q ) + ^ Im Vll 

I irv-j p L 


(B66) 


Collecting terms and evaluating the imaginary parts of the equation results 
in the pressure coefficient given in physical coordinates due to the first 
order solution given as 


= — 1 2 - ^~ Y 4=_ + — - yj x-ytanA^ 

* 1 B^ 2 y/x-ytanA 0 2 J 
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where 0 £ - ^B^+tan^A 
A = 0 H (l-ikx c ) 


* - C-x £$ 
y * n-y s 


x,y = physical coordinates 


The pressure coefficient due to the second order term may be obtained by 
using equation (B55) , equation (B63) , equation (B65) and applying the 
above procedure to yield 
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The sum of equations (B67) and(B68) will provide proper change in 
boundary conditions at the corner of a leading edge control surface. These 
expressions need to be modified such that the planform boundary conditions 
are satisfied at the edges of the planform while maintaining the strength of 
the singularities along the side edges of the control surface. The functions 
that satisfy these constraints are described in the main section of the 
document following the loading function developments of Landahl. 
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APPENDIX C 


A CAUTION REGARDING PLANFORMS WITH DISCONTINUOUS EDGES, AND A PROVISION 
FOR INCLUDING EFFECTS OF AIRFOIL THICKNESS VIA LOCAL LINEARIZATION 

Suggestions are made for modifying analytical procedures used in obtaining 
loadings on discontinuous planform shapes, and for modifying lifting surface 
boundary conditions w/V such that the physical flow pressure distributions 
are more readily simulated by the theoretical pressure distributions. 

Modification of Planforms Having Discontinuous Shapes* 

Analytical difficulties may be encountered in predicting the aerodynamic 
loadings on planforms having analytic discontinuities in the planform shape 
definition. Figure 50 represents a typical analysis planform where first 
derivative discontinuities in the shape definition exist and may provide 
inconsistent results if the spanwise distribution of collocation stations 
contain a control point located at or very near the spanwise discontinuity 
station. 



Figure 50.- Typical Analysis Configuration with a Discontinuous Shape Definition 


The mathematically derived downwash sheets obtained from an evaluation of 
the downwash Integral will exhibit singularities whenever collocation 
stations are located at the planform spanwise discontinuity stations. The 
singularities are due to the over-simplified manner in which the pressure 
distributions are assumed to exist over the planform. The assumed 
pressure functions are distributed continuously over the surface in such a 
manner that the planform is completely covered regardless of the analytic 
continuity of the planform definition. If the planform definition con- 
tains sharp breaks or discontinuities in its first derivative definition, 
then the pressure distributions and resulting spanwise loadings will also 
exhibit discontinuities at the same spanwise stations. It is the discon- 
tinuous first derivative of spanwise loadings that cause the numerical 
problems and produce singularities in the predicted downwash sheets. The 
spanwise loading discontinuities may be removed by using a suitable pres- 
sure distribution that is analytically continuous in its spanwise direction. 
However, no such distributions have been formulated nor used within the 
present program development since it is obvious that there are other means 
available for obtaining reasonable results other than overemphasizing the 
effects of localized flow in regions of planform discontinuities. The 
problem of downwash singularities may be circumvented by smoothing the 
planform discontinuities or by keeping the spanwise location of collocation 
stations well away from the troublesome regions. However, numerical prob- 
lems do exist and the analyst should be cognizant of the potential problem 
and to what extent that the downwash distributions may be distorted in the 
regions of planform discontinuities. 

The extent of downwash distortion may be determined by evaluating the down- 
wash integral equation for the steady flow case given the following form: 


W(x, y) 



In this example, f(n) is taken to be the spanwise pressure function 
corresponding to an ellipse, f(n) =)Js 2 - n 2 ; and g(s,n) is the 
chordwise pressure function given as 

s(C.n) • [ x t (n) - 5 ]' /2 / U - x t (n) ] 1/2 


(Cl) 
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Also, the plapform being evaluated is a delta wing having a discontinuous 
planform definition at the center line. 


Performing the indicated differentiation with respect to TV in equation (Cl) 
and retaining only the essential terms that contribute to the singularity 
problem results in the following expression: 


ff(x, y) = 


1 1 r «•> “h-t; 


x-C 


-S) 2+ B 2 (y-n) 2 


]«) 


dn 

y-n 


(C2) 


The term has a finite discontinuity at the planform center line 

an 

having equal but opposite signs on either side of the center line. 


To shorten the expression let the chordwise integral be defined as: 

H(n) " /* = ] d 5 (C3) 

\ /(x-e) 2 +e 2 (y-n ) 2 

And let H a (n) be H(n) when the left hand side derivative of 3 ^ - ,n ^ 
used in equation (C3) ; also let H r (n) be H(n) of equation (C3) using 

the right hand side derivative of 

Then the shortened form of equation (C2) becomes: 
s 
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Equation (C4) may be evaluated by adding and subtracting the function 

H # (n) over the integration interval o< n ^s. 
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The second and third integrals of (C5) are evaluated by making a con- 
tinuous extention of the definition of ^ across the centerline 
to the right hand tip atn» s. 

The integrals may be combined to provide the expression 


- f H^n) £ - [ 

-s J o 


j £ * f 

•'A •'a 


V") 


Singularities at n = y are removed from the integral and evaluated 
separately as shown by the following: 
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The only terms containing singularities (with the restriction that the 
downwash stations are contained in the interval -s^y*? s) are the terms 
having an integral form of 

s 

f d n 
J y - n 
o 

Equation (C7) may then be expressed as having singular and non-singular 
terms that take the form of: 


W(x,y) - [ H p (y) - H £ (y) ] ln|y| + Regular Terms 
and the limiting value of W(x,y) tends to infinity as y— ►o. 
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Thus, the downwash sheet will exhibit unrealistically large values of 
downwashes whenever a downwash chord is placed near a planform disconti- 
nuity station and is due solely to the manner In which the pressure load- 
ings are assumed to exist on the boundaries of the planform. 


There Is an obvious solution to this problem (see reference 8) in that 
the derivative of the spanwise loadings may be made continuous by smoothing 
the planform discontinuities. An example of how smoothing of planform dis- 
continuities affects the downwash distribution is shown in figure 51. 


T ^ 

SECOND DERIVATIVE 




The clipped delta planforro shown In figure 51 has been modified in the 
region of the apex at the centerline to produce a continuous first deriva- 
tive definition of the planform across the centerline. The method se- 
lected to smooth the shape definition is one of applying a function 
having a continuous first derivative that is zero at the centerline and 
is also continuous in the second derivative at the spanwise matching sta- 
tion. This particular smoothing function provides derivative continuity 
with the least distortion in planform definition. Downwash values ob- 
tained at spanwise stations along the 50% chordline are shown in the 
lower part of figure 5i . The downwash distribution near the planform 
centerline does become singular for the pointed apex case. However, the 
downwashes take on more realistic values whenever the apex discontinuity 
is removed by applying a small amount of smoothing. 

The spanwise distortion in downwash sheet is only slightly affected by 
redefining the planform shape. However, caution should be exercised in 
applying large amounts of smoothing since the redefined shape may provide 
a poor approximation of the original planform shape. 

The above discussion has been presented to identify a potential problem 
that may arise in predicting the loadings on discontinuous planforms and 
to suggest a method to alleviate the problem. The suggested smoothing 
technique is not available within the present computer program, however, 
an alternate method is suggested such that reasonable results may be ob- 
tained without redefining the planform shape. 

An alternate method suggested by other authors is one of applying engineer- 
ing judgement to the placement of spanwise collocation stations in regions 
of planform discontinuities. The downwash distribution shown in figure 51 
indicates that the distributions are only slightly affected at small dis- 
tances away from the discontinuity station. The amount of distortion will 
be the greatest near the centerline of a highly swept planform since this is 
the location where the largest difference in the spanwise loading derivative 
exists. Planforms having leading and trailing edges defined by a combination 
of straight line segments for which differences in spanwise slopes are small 
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will produce even smaller distortions in the downwash description than Is 
observed for the centerlines station. Consequently, It Is suggested that 
the placement of the spanwise collocation stations (downwash chords) is 
such that the stations are located at least a small distance away from any 
planform discontinuity station. 

Although application of planform smoothing or positioning of collocation 
stations may appear to be an artificial method to bypass the source of the 
problem, these methods will provide reasonable results for the aerodynamic 
loadings without being overly encumbered by localized flow conditions that 
conti bute only in a small way to final load definitions. 

Suggested Modification of Boundary Conditions w/V 

From an operational standpoint, it appears that if a local linearization Is 
used, that is if the linearized boundary conditions w/V are modified to Include 
local velocities due to airfoil thickness the resulting theoretical pressure 
distributions will simulate the physical flow conditions more accurately. 
Lifting surface theory solutions using assumed mode-kernel function approach 
usually provide results that correlate well with experimental results. An 
example of the correlation of theoretical -experimental results is presented 
for the analysis planform shown in figure 52. 



Figure 52.- Experimental Planform of NACA RM A51G31 (dimensions in cm) 
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Chordwise pressure distributions were obtained at seven spanwise locations / 
as indicated in the figure. The comparisons of theoretical -experimental j 
results are shown in figure 53 for the pressure chord located at y/s = .?31 
for M = 0.80 and a = 4°. 



Figure 53.- Experimental and Zero Thickness Theoretical Pressure Distribution 
at y/s = 0.831 , M = 0.80, a - 4° 


The resulting theoretical pressure distribution was obtained using the 
standard linearized boundary conditions applicable for zero thickness lift- 
ing surfaces. It appears that a reasonable correlation is obtained on the 
average in that the theoretical distribution does approximate the experi- 
mental values in an average sense over the length of the streamwise chord. 
Theoretical pressures forward of the midchord are smaller than experimental 
values and are greater than experimental values for stations aft of the 
midchord. 
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The analysis was revised using a modification on the lifting surface boundary 
conditions such that the local velocity due to symmetrical thickness distri- 
bution is used instead of applying the uniform velocity distribution of a 
zero thickness lifting surface. 


That Is, the boundary condition of the integral equation^was obtained from 
the definition 


fi (x,y,o) =}§| = }[at i+ 3x i ||] 


where z. represents the displacement of the surface and is a function of time(t). 

J 


Ifi 

3X 


the slope of the surface at (x,y). 


3 X 

H* is usually taken to be equal to the remote velocity V, and for a zero 
thickness airfoil section this may be correct. However, physical experiments 
are conducted on finite thickness airfoil sections having local velocities 
that are not uniform over the chord length. 


Consequently, the ^ term should represent the local flow velocity ^i oca i at 
collocation stations in making a comparison between theoretical and real 
flow results. 


The modified boundary conditions applicable to finite thickness airfoil sec 
tions are then defined as: 


W (x,y,o) 

3 


= j. r Ifi 

V L 3t 


dz • 3Z. 


ax Local 


] 


\ocal 1S defined as being the steady streamwise velocity distribution that 
differs from V due to thickness effects only, and may be obtained from ex- 
perimental results or by using the theoretical distributions of reference 17. 


Theoretical results shown in figure 54 represent analyses using zero thick- 
ness and finite thickness boundary conditions and it is evident that in 
this case more accurate simulation of the physical flow results are obtained 
using local boundary conditions that incorporate finite thickness effects. 
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Theoretical pressure distributions forward of the midchord station approach 
the experimental values in a better fashion than do the results using zero 
thickness boundary conditions. Also* the theoretical distributions near 
the trailing edge predict the change in sign in pressures as Indicated in 
the experimental results. 



x/C 

Figure 54.- Comparison of Pressure Distributions Predicted by Zero Thickness 
and Finite Thickness Boundary Conditions at a Semi-Span Station 
y/s = 0.831, for M = 0.80 and a = 4° 


Although no large changes appear In the distribution using the modified 
boundary conditions, it does appear that the physical flow conditions are more 
accurately simulated at least for this steady flow test configuration. 

Validity of applying local linearization of boundary conditions within the 
unsteady case may be questioned on the basis that the kernel function has not 
been changed to reflect the local linearization effects. 
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However, a numerical experiment has been accomplished to evaluate the 
effects of- local linearization on the unsteady pressure distributions using the 
test configuration shown In figure 55. 


0 . 60 m 



Figure 55 - Experimental Planform of Slde-by-SIde Control Surface Arrangement 
Showing Spanwlse Location of Pressure Measuring Stations. 


Results of this investigation are shown in figure 56 for the in-phase part 
of the resulting pressure distributions and the out-of-phase part is presented 
In figure 57. 

The theoretical pressure distributions of the In-phase part approach the 
experimental values in a better fashion than do the results using zero thickness 
boundary conditions. The out-of-phase theoretical distributions remain almost 
Invariant over the control surface and are changed only by very small amounts 
forward of the hingeline. 
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Figure 56 - In-Phase Pressure Distribution for Both Flaps Oscillating with 
the Same Phase and Amplitude. 6 f “ 0.82*,k« 0.372 , M ■ 0.0 



Figure 57 - Out-of-Phase Pressure Distribution for Both Flaps Oscillating 
with the Same Phase and Amplitude. 0.82 » k* 0.372, M" 0.0 




It should be noted that the validity of including local linearization 
effects within the boundary conditions can not be generalized on the basis of 
this single example. 

However, it appears that the physical flow conditions may be more 
accurately simulated and may contribute significantly to the successful design 
of energy absorbing Stability Augmentation Systems provided that good correla- 
tion between experimental and theoretical results can be established over the 
frequency range of interest. 

A program option is available for modifying the boundary conditions 
within the analysis and the data format along with limitations are given in 
reference 4. 
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